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Abstract: Infrared power corrections for Minkowskian QCD observables are ana- 
lyzed in the framework of renormalon resummation, motivated by analogy with the 
skeleton expansion in QED and the BLM approach. Performing the "massive gluon" 
^ ■ renormalon integral a renormalization scheme invariant result is obtained. Various 

regularizations of the integral are studied. In particular, we compare the infrared 
^ ■ cutoff regularization with the standard principal value Borel sum and show that they 

■ yield equivalent results once power terms are included. As an example the average 

thrust (T) in e~^e~ annihilation is analyzed. We find that a major part of the dis- 
crepancy between the known next-to-leading order calculation and experiment can 
be explained by resummation of higher order perturbative terms. This fact does 
not preclude the infrared finite coupling interpretation with a substantial 1/Q power 
term. Fitting the regularized perturbative sum plus a 1/Q term to experimental 
data yields ^^(Mz) = 0.110 ±0.002. 
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1. Introduction 

Power corrections in QCD have been the subject of many interesting theoretical 
developments in the recent years [1], especially concerning observables that do not 
admit an operator product expansion. Of particular interest in this respect are 
event shape variables in e+e~ annihilation. For these observables the available [2] 
next-to-leading order perturbative calculations in a standard choice of renormal- 
ization scheme and scale are found not to agree with experimental data [3]. The 
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discrepancy was originally bridged over by Monte-Carlo simulations which account 
for non-perturbative "hadronization corrections". The data can also be fitted by 
the next-to-leading order perturbative result plus a power correction that typically 
(but not always) falls as 1/Q, where Q = ^/s is the center of mass energy. Nev- 
ertheless, this situation is not satisfactory, especially because the non-perturbative 
correction, which is not under control theoretically, is numerically quite significant. 
Consequently, much theoretical effort has been invested in the last five years in un- 
derstanding the source of these power corrections. This effort turned out to be quite 
fruitful [4]-[ll]: a "renormalon phenomenology" has been developed [1], where infor- 
mation contained in perturbation theory is used to determine the form of the power 
correction, while its normalization is left as a non-perturbative parameter, which is 
determined by fitting experimental data. 

On the other hand, the very same observables appear to have significant sub- 
leading perturbative corrections and large renormalization scheme and scale depen- 
dence when calculated up to the next-to-leading order, as discussed in ref. [12, 13, 3]. 
This observation raises concern about the rehabihty of the results obtained in the 
usual procedure, where a perturbative expansion truncated at the next-to-leading 
order is used^ as a starting point for the experimental fit [3]. Unfortunately, the full 
next-to-next-to-leading order calculation for these observables is not yet available. 

In this paper, we adopt the view that the most important corrections are related 
to the running of the coupling, in the spirit of the BLM approach [15]. Following 
[16], we assume that the perturbative series can be reshuffled in a (yet hypothetical) 
"dressed skeleton expansion" built in analogy with the Abehan theory, where each 
term is by itself renormalization scheme invariant. This expansion was the original 
motivation behind the BLM approach (see ref. [17] for further discussion). In the 
Abehan theory, the skeleton expansion is well defined. The first term corresponds to 
an exchange of a single photon, dressed by all possible vacuum polarization insertions 
which build up the Cell-Mann Low effective-charge; the second term corresponds to 
an exchange of two dressed photons and so on. The expansion coincides with the 
standard expansion in a for a conformal theory where the coupling does not run. 
The leading skeleton term can be written compactly as a "renormalon integral" - 
an integral over all scales of the Gell-Mann Low effective-charge times an observable 
dependent function which arises from the one- loop Feynman integrand and represents 
the momentum distribution of the exchanged dressed photon. When the leading 
skeleton term is expanded in terms of the coupling up to large enough order, both 

^For event shape distributions close to the 2-jet hinit, large Sudakov logarithms related to the 
emission of soft and coUinear gluons appear explicitly in the perturbative coefficients. Due to 
multiple emission such contributions persist at high orders and make the fixed order calculation 
useless. In ref. [14] it was shown how these large logs can be systematically resummed to all orders. 
In this work we concentrate on average event shape observables, where these contributions are 
presumably not important. 



2 



large and small momentum regions give rise to factorially increasing perturbative 
coefficients. These are the renormalons [18]-[21] which are believed to dominate the 
diverging large order behavior of the full perturbative series. This way renormalons 
carry information about the inconsistency of perturbation theory and thus on possible 
non-perturbative power corrections arising from the strong coupling regime. 

Having no straightforward diagrammatic interpretation in the non-Abelian case 
the "skeleton expansion" is still just a conjecture. There is an obvious difficulty 
due to gluon self-interaction vertices. This difficulty can be hopefully resolved by 
separating contributions from such diagrams into several different skeleton terms in 
a gauge invariant manner. Ref. [22] gives a concrete suggestion for a diagrammatic 
construction of the skeleton expansion in QCD, based on the pinch technique. How- 
ever, more work is required to establish it beyond the one-loop level. On the other 
hand, the structure of the single gluon exchange term ( "leading skeleton" ) is strongly 
motivated by the large Nf limit where the Abelian correspondence is transparent. 
The non-Abelian analogue of the renormalon integral, obtained through the "Naive 
Non-Abelianization" [23] procedure, was used extensively in the last few years in 
all order perturbative resummation (e.g. [24]-[27]) and parametrization of power 
corrections (e.g. [6]-[ll] and [16]). 

In the present work, we perform an analysis of infrared power corrections to- 
gether with all order perturbative resummation of renormalon-type diagrams in the 

"massive gluon" approach [24, 25, 8]. We show that these two aspects of improving 
the standard perturbative calculation cannot be dissociated and must be performed 
together. 

We study one specific example: the average thrust in e~^e~ annihilation, where 
there is a priori evidence for both large perturbative corrections and strong 1/Q power 
corrections. We perform renormalon resummation at the level of a single gluon emis- 
sion, emphasizing the renormalization group invariance of this procedure. We discuss 
in detail the ambiguity between the perturbative sum and the non-perturbative in- 
frared power corrections. We also address the complications that arise when applying 
the inclusive "massive gluon" resummation to not-completely-inclusive Minkowskian 
observables, such as the example at hand. 

The paper is organized as follows: in sec. 2 we describe the specific assumptions 
we make concerning the "dressed skeleton expansion" in QCD and the immediate 
consequences that follow. We also compare the "skeleton expansion" with the BLM 
scale fixing procedure. In sec. 3 we discuss various natural regularizations of the per- 
turbative sum, concentrating on Minkowskian quantities. In sec. 4 we show that the 
regularization independence of the full QCD result is achieved only by adding to the 
regularized perturbative sum explicit power terms. We also discuss the connection 
with the infrared finite coupling approach [7, 8] . In sec. 5 we present the application 
of the method to the average thrust. Sec. 6 contains our conclusions. 
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2. The "dressed skeleton expansion" and BLM 



Consider first a generic Euclidean quantity D{Q^), which has the perturbative ex- 
pansion 

^pt(<5^) = roa + rio^ + • • • (2.1) 

where a = ^ is the couphng at scale nj^ in (say) the MS scheme. Since the leading 
coefficient Tq is flavor independent and the next-to-leading coefficient ri is linear in 
Nf, one can write 

-Po I log — 



where 



n = ro 



Po 



71 



5i 



(2.2) 



1 /II 



is the one-loop {3 function coefficient, /3(a) = da/d\n fj.'j^ = —Poo^ + • • • , and 71 and 
Si are flavor independent. We shall assume, in analogy with QED, that DpT.{Q'^) has 
the "dressed skeleton expansion" 



D,,{Q') = D^,''{Q') + Dr{Q') + 



(2.4) 



where D^'^ is the contribution of a single dressed gluon exchange, comes from a 
double exchange, etc. This means in particular that we assume for D^'^ the repre- 
sentation 

d,^(e)^n{kyQ') (2.5) 



A:2 



where $_d is the "momentum distribution function" [26] which depends on the ob- 
servable under consideration {D), whereas would be given by 

DriQ') ^f^^ «pt(^?) -a^kl) ^Ukl/Q', (2.6) 

and so on. The physical "skeleton coupling" appearing in both (2.5) and (2.6) 
does not depend on the observable under consideration. It is supposed to be uniquely 
determined-'- and is a priori different from the renormalization scheme coupling a = 
a(/i|.). We can therefore consider its (a priori non-trivial) expansion in the renor- 
malization scheme coupling 



-Po (log 

V f^R 



Ci 



di 



Using eq. (2.7) in eq. (2.5) we obtain 



(2.7) 



(2.8) 



■t-In QED, it is the Gell-Mann Low effective charge. 
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with To = 00 and 



r? = 0o 



-/?o log ^ + -i + ci + (ii 



(2.9) 



where (pi are the log-moments of the momentum distribution function, 

The only source of Nf dependence in the next-to-leading order coefficient ri are 
vacuum polarization corrections which dress the exchanged gluon, namely corrections 
that are fully included in Dq'^{Q'^). It therefore follows that the "Abehan" parts 
proportional to /3o in ri and are the same, i.e. 

J^oTi = 0i + 0oCi. (2.11) 

The next observation is that the "leading skeleton" term D^'^{Q'^) is a renor- 
malization scheme invariant quantity, just as the "skeleton coupling" apT(/c^) itself. 
The question whether the approximation of DpT^^Q"^) by L'q^(Q^) is a good one thus 
has a renormalization scheme invariant meaning. In particular the leading order 
correction, which is 0{a^), can be written as 

Dp,(g2) ~ D^-^Q^) + - r°) + • • • (2.12) 

where ri — is renormalization scheme invariant. Note that 

ri - r? = ro {S, - d,) (2.13) 

and the r.h.s. can be identified as the difference between the next-to-leading coeffi- 
cients which remain after performing BLM scale setting [15] in the two renormaliza- 
tion group invariant quantities DpT,(Q^) and Dq'^^Q^) respectively. Such a difference 
is known [15, 28] to be renormalization scheme invariant, although 5i and di are 
separately scheme dependent. 

It is interesting to compare the "skeleton expansion" approach which we take 
here with the standard BLM scale setting procedure [15]: 

1. Eq. (2.11) is equivalent to the statement that the BLM scale for the quantities 
Dpt{Q'^) and Dq'^{Q'^) is the same: 

/^Lm = exp + ci^ = exp(7i) (2.14) 

with (f)i defined in (2.10). Note in particular that if BLM is apphed in the 
"skeleton scheme" , where = and then the BLM scale is the center 
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[26, 29] of the momentum distribution^ ^oik'^/Q^), 

/^BLM = Q exp — 




(2.15) 



i.e. it is the average virtuahty of the exchanged dressed gluon. Using the 
BLM scale, the leading order term roadilj^^) can be viewed [26, 29] as an 
approximation to the entire all-order sum Dq'^{Q^). This approximation is 
good if ^D{k^/Q^) is narrow, and it is exact for 



2. The usual justification of BLM (in a generic scheme) relies on the assumption 

that Si is small and so the large /5o piece — /5o ^log ^ + 71^ dominates the full 

next-to-leading coefficient ri in eq. (2.2). This depends on the renormalization 
scheme, scale and Nf. On the other hand, in the "skeleton expansion" approach 
there is no scale setting to perform, and so the accuracy of approximating 
-Dpt(Q^) by Dq'^{Q^) is controlled by the magnitude of a scheme invariant 
coefficient (2.13); the issue now is whether di is a good approximation to 5i. 
Note also that if the quantities in eq. (2.13) are computed in the "skeleton 
scheme" then' di — and the scheme invariant parameter that controls the 
accuracy of the "skeleton expansion" can be identified as the standard BLM 
coefficient ro Si. 

As mentioned in the introduction, it is not yet clear whether a "skeleton expan- 
sion" exists in QCD. Thus we do not know the identity of the "physical skeleton 
couphng" OpT- We do know, however, that in the Abelian hmit the "skeleton cou- 
pling" should coincide with the V-scheme and so the Abelian coefficient ci is deter- 
mined. For instance, in the MS scheme, Ci = — |. In order to perform the "leading 
skeleton" rcsummation in practice wc need to specify also the non-Abelian coefficient 
di. We shall therefore consider in this paper three schemes which share the same Ci 
mentioned above and differ by di (all quoted values arc in MS ) : 

a) The "gluon bremsstrahlung" coupling [30], where d\"'"^ = 1 — This cou- 
pling was used in [7, 8] for parametrization of power corrections to event-shape 



b) The skeleton effective charge found in [22] using the pinch technique, where 



§In order to interpret <!>£) as a distribution function it has to be positive definite. Although no 
general argument is currently known, it turns out that in practice is positive definite for almost 
all investigated physical quantities [29] . 




(2.16) 



variables. 



= 1. 



^In QED di is always zero. 
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c) The V-scheme coupling [15] , defined by the static heavy quark potential, where 
dX = -2. 

One might worry that without a precise identification of the "skeleton coupling" we 
introduce back some kind of renormalization scheme dependence. We shall see that 
in practice (sec. 5), the inclusion of a next-to-leading order correction as in (2.12) 
effectively compensates to a large extent for the ambiguity in di. 

3. Regularization 

Wc begin with the observation that the integrals representing the terms DJ'^ in the 
"skeleton expansion", e.g. Dq'^ of eq. (2.5), are ill-defined since they run over the 
Landau singularity. This is the way the infrared renormalon ambiguity of the pertur- 
bative sum appears in the framework of the "skeleton expansion" . It is possible to 
make these renormalon integrals mathematically well-defined by specifying a formal 
procedure to avoid the Landau singularity. However, this would not cure the physical 
problem of infrared renormalons: some additional information associated with large 
distances is required in order to obtain the full QCD result from the perturbative 
one. 

Still, in order to use the "skeleton expansion" in practice one is bound to define 
Df"^ somehow. We shall concentrate in this paper on the single dressed gluon term 
Dq'^ and refer to the regularized integral of eq. (2.5) as the "sum of perturbation the- 
ory" . We shall consider various possible regularization prescriptions: Principal Value 
(PV) Borel summation. Analytic Perturbation Theory (APT) approach ( "gluon mass 
integral"), infrared cutoff method, and truncation of the perturbative series at the 
minimal term. Any two such prescriptions differ just by power terms. Thus, the 
"dressed skeleton" approach leads us implicitly to the subject of power corrections. 

It is clear that physics does not depend upon the definition used. In practice, hav- 
ing no way to calculate the non-perturbative contribution we shall just parametrize it 
properly and fit the data by the "sum of perturbation theory" plus a power term. The 
prescription dependence of the two separate contributions cancels by construction. 
While in principle any regularization can be used, some may be more illuminating 
then others. As we shall see in sec. 3.3 and later in sec. 4 the infrared cutoff regu- 
larization is of special interest, allowing to separate at once short vs. long distance 
physics and perturbative vs. non-perturbative physics. 

The purpose of this section is to derive the general relations between different reg- 
ularization prescriptions which are useful in the analysis of Minkowskian observables. 
We begin by discussing the application of the "skeleton expansion" to Minkowskian 
observables concentrating on the "leading skeleton" term. We review the Analytic 
Perturbation Theory (APT) or "gluon mass" [24, 25, 8, 16] integral which seems to 
be the most convenient regularization for a practical calculation of the perturbative 
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sum (sec. 3.1). We then (sec. 3.2) discuss the power corrections that distinguish be- 
tween the APT integral and the principal value Borel sum [24, 25, 16]. In sec. 3.3 we 
explain how an Euclidean infrared cutoff can be introduced for Minkowskian observ- 
ables and derive the relation between the cutoff regularized perturbative sum and 
the APT integral. In sec. 3.4 we generalize the explicit formulae of sec. 3.2 and 3.3 to 
the case of a two-loop "skeleton coupling", and finally, in sec. 3.5 we comment of the 
comparison between different regularizations. Appendix A describes an alternative 
computation method for the cutoff regularization in terms of the APT integral. 

3.1 Minkowskian quantities and the APT integral 

We assume that a "skeleton expansion" such as (2.4) can be constructed also for 
Minkowskian observables. For a Minkowskian observable which is related by a dis- 
persion relation to an Euclidean quantity, linearity of the dispersion relation implies 
that the expansion will take the form 

R.t{Q') = R7{Q') + R7iQ') + ■■■ (3.1) 

where Rq'^ is the leading "dressed skeleton" which is related to Dg"^ by a dispersion 
relation [16]. Similarly, should be related by a dispersion relation to of eq. 
(2.6), and so on. If there is no dispersion relation with an Euclidean quantity, the 
existence of a "skeleton expansion" is more doubtful. In particular, as we discuss 
in sec. 5, there is presumably no way to replace the entire perturbative series of 
not-completely-inclusive observables such as weighted cross sections by a "skeleton 
expansion" . 

Let us concentrate now on the "single dressed gluon approximation" Rq"^. It 
turns out that Rq'^ cannot be expressed in the "Euclidean" representation of eq. (2.5) 
as an integral over the space-like coupling ap^. Instead, R^^ has a "Minkowskian" 
representation'! in terms of the time-like discontinuity of the coupling [24, 25, 26, 8], 

Rq^ = Rapt with 



- r^c^(/^^)^«(/^vQ^) (3-2) 

Jo jJ.'^ 

where the reason for the name APT shall be explained below and 

PPT(At^) = ^^^^^ |apT(-A*^)} = ^ Opt {-jJi^ + ie) - (-^^^ - ie) (3.3) 

"This representation applies only to inclusive enough quantities which do not resolve the de- 
cay products of an emitted gluon [24, 25, 8]. Then the time-like "skeleton coupling" a^^ can be 
reconstructed from the higher order terms related to the gluon decay. 
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is the time-like "spectral density". The corresponding Minkowskian "effective cou- 
pling" (/x^) is defined by 

=Pp.(/.^). (3.4) 



dln/i^ 

Note that this time-like coupling is well defined and finite at any /i already at the 
perturbative level, as opposed to the space-like coupling that in general has Landau 
singularities. For instance, in the one-loop case 

apT(A;^) = (3.5) 
Po log 

has a Landau singularity while the corresponding time-like coupling. 



1 1 

arctan 

2 TT 




(3.6) 



is finite for < /x^ < oo and has an infrared fixed-point: (0) = l//3o- 

The "characteristic function" J^r in eq. (3.2) is computed from the one-loop 
Feynman diagrams with a finite gluon mass /i, and = —dJ-'u/ din /i'^ . It is usually 
composed of two distinct pieces 

where is the sum of a real and a virtual contribution, while J-'^^ contains 

only the virtual contribution, and may vanish identically as in the case of thrust. 
This property prevents [26, 16] a representation of i?Q^ similar to eq. (2.5) to be 
reconstructed from eq. (3.2) using analyticity. 

In the Euclidean case, the characteristic function is made instead of a single 
piece and satisfies a dispersion relation. Its discontinuity is [26, 31] nothing but the 
Euchdean "momentum distribution function" $d of eq. (2.5) 



2m 





(3.8) 



However, 



D^,M) ^ a-(/.2) j^^ii^yq^) (3.9) 

Jo 



dfi^ 
J' 

differs from D^'^ of (2.5) by power terms: while Dq'^ is ambiguous involving inte- 
gration over the Landau singularity, -Dapt is well defined thanks to the non-trivial 
infrared fixed point in a^j^(/x^). e.g. (3.6) in the one-loop case. 

One can use analytic continuation to express /^apt in the Euclidean representa- 
tion. This yields 

D,MQ') J-^ «apt(A;') ^Dik'/Q') (3.10) 
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where defined through the following dispersion relation 



2 



T-T^KIil^') (3.11) 



/o (/i^ + P)2 

which implies in particular the absence of Landau singularity. The name Analytic 
Perturbation Theory [32] is attached to this coupling since it has, by definition, a 
physical analyticity structure. The APT coupling differs from the usual perturbative 
coupling by power terms 

OaptI^^) = apT(/c^) + (5aAPT(^^) (3-12) 

with 

5a^,.{k^) = Y.i-^T^'K[^] (3.13) 



n=l 

and in the one-loop case 



5a^..{k^) = --^^A-r- (3.14) 

PO ;p - 1 

Returning now to Minkowskian quantities, we keep the subscript "APT" in 
eq. (3.2) in order to stress the fact that this "gluon mass integral" is well defined and 
is different from the corresponding Borel sum R^'^ . Rapt and Rq^ share by definition 
the same perturbative expansion. The two differ by power corrections which are in 

general ambiguous, and thus i?APT can also be looked at as a particular regularization 

•PI 





of Rr 



3.2 APT vs. Borel sum 

It is useful to derive an expression of the Borel sum** i?pT in terms of Rapt, following 
[24, 25]. The "renormalization scheme invariant" Borel representation [33, 34, 35] of 
eq. (3.2) can be written as [16] 

RMQ') = dz -a^^z) ^ TniiJ^/Q^) exp l^-z(3^ In (3.15) 

where d^siz) is the (renormalization scheme invariant) Borel image of the Minkowskian 
coupling 

KI{l^^)=j^ (i^exp ( -^/3oln — j a^s{z). (3.16) 
a^s{z) is related to the Borel image a{z) of the space-like coupling 

/■oo / /c^ \ 

apT{k'^) — J dzexp l—zPoln — j d{z) (3.17) 

**The "leading skeleton" subscript is dropped in the following for simplicity: RpT refers from 
now on to the Borel sum of the "leading skeleton" (-Rg^ in the previous section). 
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by the relation [33] 



(3.18) 



We define A to be the Landau singularity of and assume that the Borel integrals 
converge for I? and /x^ larger then and that the Borel representation of i?pT 
converges for large enough Q^. Eq. (3.15) is obtained by substituting eq. (3.16) into 
eq. (3.2), and interchanging the order of integrations. 
Next we choose a separation scale > A and write 



i?p,(g2) = i?7(g^) + i?^-(g2) 



where 



POO r 

Jo Jo 



2 ^ 

A2 



(3.19) 



(3.20) 



contains the infrared renormalon ambiguities, and 



POO r 

Rl^iQ') = / dz a,,{z) / 

Jo Ju 



2 Mi^VQ^) exp -^/3oln 



A2 



(3.21) 



is unambiguous. Since /ij > A, we can replace a^J" (/x^) by its Borel representation 
for II > II I, and therefore 



Jlxf II 



where 



Jo 



11^ 



We thus end up with 



where [16] 



RpriQ^) — -RAPT(g^) — 5-RAPT(g^) 



SR,,,{Q') = Rf-{Q')-Rl-{Q'). 



At the one-loop level 5-Rapt can be expressed quite elegantly as [24, 25] 

1 

'Wo 



SR 



APTli.loop 



(3.22) 
(3.23) 

(3.24) 
(3.25) 

(3.26) 



The principal value Borel sum -Rpt|pv is then obtained by replacing 5-Rapt by its real 
part in eq. (3.24). The imaginary part of SR^pt provides a measure of the summation 
ambiguity. 

Note that although R^'^ and R^ in (3.25) are separately iij dependent, their 
difference SRapt does not depend on This formula is of practical utility when 
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only the small j-i^ /Q"^ expansion of J^r is known analytically (e.g. the average thrust; 
see sec. 5): -Rapt can then be obtained from (3.2) by numerical integration over JFr, 
whereas (5i?APT, which involves only scales below /x/, can be evaluated at large Q"^ 
using the small fJ^ /Q^ expansion of Tr. 

Consider indeed the contribution of a generic term in the small [i^ jQ^ expansion 
of Tr^ namely assume 



.F„(„W)-J-h(0)=^ \Bt'\oi^ + Ct'\ (3.27) 



where n is not necessarily integer. Eq. (3.27) implies 

j=K{^'m = -n ^ (l^r log ^ + 4"^ - ^ j ■ (3.28) 

We then obtain from eq. (3.20) 

VQV ^ "'^ ' (i- 

where Zn= and from eq. (3.23) 

R <JXQ') = 



^ dz ~a^,iz) -2 exp -zPo In % (3.29) 

n Jo M_zr V A^/ 



n^^F^ log^C^V). (3.30) 



In both (3.29) and (3.30) the combination in front of the first integral (the single 
Borel pole integral in (3.29)) is the same as in J^r. The double Borel pole integral 
in (3.29) originates from the logarithmic term in (3.27) and it can be obtained by 
differentiating the single pole integral with respect to n. 
One deduces the generic contribution to SRapt at large 



/A2 



2 



n n 



(3.31) 



with 



^n^l ^^[j-2j - i ^-(^) (3.32) 
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and 

(x^)" 'V -/."""^ "-w 

where we have used the fact that SR^^'^ does not depend on ni and thus have set 
Hi = A. Note that 6^ is related to by 

The integrals (3.32) and (3.33) are in general difficult to calculate. However, 
in the one-loop case {a{z) = 1) we can obtain the final results for 6„ and 6^ by 
comparing (3.31) with (3.26): 

^'n|,,„, = ^e±^- (3.35) 

and (deriving with respect to n) 

&nli-ioop = ^e±-"(l±i7rn). (3.36) 

Eq. (3.35) can also be checked for n integer with the observation [16] that the 6„'s 
of eq. (3.32) then coincide with the coefficients in the large k'^ expansion of SaApj-ik'^) 
defined in eq. (3.13). 

3.3 Infrared cutoff regularization 

Let us now turn to the infrared cutoff regularization. Consider first an Euclidean 
quantity D^Q"^). A natural regularization of eq. (2.5) is obtained by introducing an 
infrared cutoff /// (above the Landau singularity) 

/OO flip' 
^ ^ ap,(A;^) ^nik-'m (3.37) 

which removes the long distance part of the perturbative contribution 

= Dp^{Q')-D-^{Q'). (3.38) 

The regularized perturbative sum (3.37) is fully under control in perturbation theory, 
provided /xj/A is large enough. One can then safely evaluate the "leading skeleton" 
term using a one or two loop approximation to the running "skeleton coupling" . The 
accuracy of this approximation can be estimated comparing the results obtained at 
the one and two loop levels. In this respect the cutoff regularization is of special 
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interest: other regularizations, e.g. the principal value Borel sum, involve properties 
of the coupling in the infrared < /c^ < /if. It is then hard to explain why the 
one-loop or two-loop couplings are of any relevance. 

For Minkowskian quantities, we have seen that a representation such as eq. (2.5) 
does not exist. One could imagine a regularization in the form (3.21). However, as 
opposed to a separation between large and small space-like momentum, the separa- 
tion between large and small time-like momentum (3.19) does not correspond to a 
separation between short and long distances. 

It is natural to define instead [45, 36, 16], just as for Euclidean quantities. 



where Rp-r is the Borel sum (eq. (3.15)), and 



R,, 



Jo 



(3.39) 



(3.40) 



Here ^ji{k'^/Q'^) is (minus) the discontinuity at /x^ = —k"^ < of the "small gluon 



mass" piece J^j^ ■* of the characteristic function (cf. eq. (3.8)) 



2m 




(3.41) 



The regularized perturbative sum i^^v defined in eq. (3.39) has a similar structure 
to its Euchdean analogue D^v of sq. (3.38). The only difference is that R^'^ (hke 
Rpt) cannot be written as an integral over the space-like "skeleton coupling" Cpt. 
For Euclidean quantities, the prescription of cq. (3.39, 3.40) reproduces eq. (3.37). 
In [16], it was shown explicitly using a Borel representation that R^'^ is indeed free 
of any ambiguity from infrared renormalons. We have 



RlliQ^ 



d^^ 

jJ? 



^7(/x^g2) exp -^/3oln^ 



where 



dk"^ 



'o k'^ -\- //^ 

is the "infrared regularized" characteristic function, and 



9 In //^ 



(3.42) 



(3.43) 



(3.44) 



The effect of the subtracted term in eq. (3.43) is to remove any potential non-analytic 
term in the small /x^ ("gluon mass") expansion of jF]^^(/i^, Q^) (now distinct from 
its large behavior, since there is a third scale /// involved), by introducing an 
infrared cutoff on the space-like gluon propagator momentum. This implies the 
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cancellation of infrared renormalons, which are related only to non-analytic terms 
[46, 24, 25, 8] in the small /i^ expansion of T{^'^/Q'^). The Borel representation 
eq. (3.42) is not so practical for concrete calculations, especially in cases (like the 
one of the average thrust which will be considered in sec. 5) where the characteristic 
function is not known analytically. In this section we describe an alternative method 
to compute R^, based on the "gluon mass" regularization (3.2). 
Using eq. (3.24) one can indeed write 



(3.45) 



with 



pPT pPT 



Rt 



(3.46) 



where the overall //./-dependence is that of -Rfj^. Ry^ should be free of any renormalon 
ambiguity. Since -Rapt in (3.45) is unambiguous, so is AR. We further know that 
R^'^ in (3.46) is unambiguous and thus it is clear that the ambiguities present in R^^ 
and R^ should cancel. We can therefore calculate AR as a sum of two well-defined 
contributions''^^ : 



a) {R^^ — R^) which is a Borel integral free of renormalon ambiguities. 

b) R^"^ which is the "small gluon mass" integral. 

Let us now calculate {R^^ — R^) explicitly for a generic term such as (3.27) in 
the small /x^ expansion of and check that it is indeed unambiguous. Using (3.41) 
we obtain 



which implies that 



sm nn 



7T 



(n) 



— Bj^ cos nn 



(3.47) 



B 



Smc TTTl 



(n) 
R 



COS Trn 



n 



'R 

n 



smc Trn 



(3.48) 



where sine x = 

rr. ' 



2 «pt(^ ) 



n 



k^ 



(3.49) 



''"''An alternative method to calculate Ai? based on the APT coupling is described in the appendix. 
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and 



4 {^^y^') - 



U?/ 



log — apT(/i;^ 



Replacing apT(A;^) by its Borel representation (3.17), one gets 



In (/X?/A^) 



1 f /i^^ 

(i^; a(z) ( "^-^o 



(3.50) 



(3.51) 



and 



(/^//A^) = / 5(2;) — exp \-z(3o\n 



A2 



(3.52) 



Note that both R^^^ (^l- (3-29)) and -Rfi^„ have simple and double renormalons poles 
a,t z — Zn = j^- As mentioned before, the double poles reflect the log-enhanced small 
fi^ behavior of J^r. 

As an example, in the one-loop case where d{z) = 1 the integrals in eq. (3.51) 
and (3.52) yield 

(3.53) 



n 



4|i-.oop= - -^e-"*^ Eil (-nt/) 



and 



1-loop 



n 
Wo 



nt/e-"*^Eii(-ni7) + l 



(3.54) 



where tj = In (/^l/A^) and the exponential integral function is defined in the complex 
plane by 

/OO fly: 
e— ^. (3.55) 

This function has a cut on the positive real axis, and thus and are ambiguous 
as expected from eq. (3.51) and (3.52). 

On the other hand, using eq. (3.18) and (3.29) one deduces 



(3.56) 



where 

and 



B<"' log ^ + 4"') p. („;/A^) - (|) " („;/a=) 



n 



dz d(z) — 
1 



z 



[sine Trn — sine ttPqz] exp —z^q In 



A2 



(3.57) 



Pi (/^?/A^) = 

/ d(z) 
Jo 



1 - ^ 



sine 7r/?o2; — cos vrn — 



sine t:(5qz — sine 7rn 



1 - ^ 



exp 



(3.58) 
-^/^oln^J 
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Explicit calculation of the limit of the integrands in p„ and for z — > Zn shows 
that the renormalon pole cancels. 

Note that the combination appearing in front of in (3.56) is the same as in 
J^R. Alternatively, we can rewrite eq. (3.56) such that the combination in front of Pn 
will be that of namely 




(3.59) 



B 



(n) , Q , 

^ log-^ + C^ - — 



R 



/ ,,2 \ " 



where 

A (m?/a^) 



oo \ 

^ dza{z)^^—^ 



smc Trn — cos 



sine ttPqz — sine tto 



1 - ^ 



(3.60) 

exp 1-zI3q In ^ J 



It is now straightforward to see that the expression in the square brackets vanishes 
in the limit z — Zn- 

In the one-loop case, where a{z) = 1, all the Borel integrals can be evaluated 
analytically. We obtain from (3.57) 



1 1-loop 



-^Re [i Eii {n {in - tj)) e"^^'^-*^)) + ^ 



arctan 



TT 



(3.61) 



The exponential integral function in the first term in (3.61) is defined in (3.55) and 
it is unambiguous for complex arguments. Wc identify (see (3.6)) the second term in 
(3.61) as the time-like coupling at the cutoff scale, namely a^^ (/^f). Note that this 
term is independent of n. Similarly, we obtain from (3.58) 



ft ( 

= -5— Re{(-^)Eil(r^(^7r-^,)) 

i-loop Po TT 



ti)\. (3.62) 



The second ingredient required for the calculation of Ai? in (3.46) is the "small 
gluon mass" integral which is defined in (3.23) and written explicitly in (3.30). 
In order to calculate Bf^J^ it is convenient to first integrate eq. (3.23) by parts to get 
an expression in terms of the discontinuity of the coupling 







^«(/^VQ')--^i?(o)JpPT(/x') 

(3.63) 

We take again a generic term (3.27) in J^r{ii'^ /Q^) and perform the integral analyti- 
cally in the one-loop case. The result is 



^ APT / 



<.n{Q'')= (3.64) 
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where hn = —pn + + Qn and /i^ = — + q}^ + g\. The functions g„ and g^, as 
well as Qn and gf^, will be given a simple interpretation below. In the one-loop case 



9n I i-ioop 
1 



1-loop 



1 

n 
Jo 



cos 7rn e 



—ntj 



(3.65) 



and 



1 sin Trn 



in I i-ioop 



l-loop 



,—nti 



1 



(tt sin nn + tj cos 7m)e 



Ei(nt7) 



(3.66) 



sm nn 

TT 



+ e~"*^Ei(ni/)n (t/ 



sm nn 



— cos Trn 



TT 



where Ei(a;) in (3.66) is the Cauchy principal value exponential integral function, 
defined for real arguments such that for negative x, Ei(x) = — Eii(— x) and for 
positive X, Ei{x) = Re {— Eii(— x)}, where Ei^ is defined in (3.55). 
Combining eqs. (3.56) and (3.64) according to 



AR = {R^;^ - i?^^) + Rf^ 
we find that the generic large contribution Ai?„ is given by 



(3.67) 



ARniQ' 



R 



2\n j^(n) 



R 



n 



(3.68) 



Finally, we can rewrite the result as indicated in the first line of (3.46) separating 
AR into i?fj^ and SR^pt = R^^ — R^- As explained above, R^^ and 6Rapt are both 
imaginary and ambiguous. Nevertheless, knowing that the imaginary parts cancel in 
AR, we can write 

ARn = Re {i?f,^4 + Re {5Rf-} . (3.69) 

This separation of AR makes a clear connection between the cutoff regularization 
and the Borel sum principal value regularization i?pT|pv, 



-^uv — -Rapt — Re {dRj^pr^} — Re {-Rfa } — -RpT|pv — R-e {-Rfft } 
We now make the following crucial identification: 

Re {-^rH^n} = 



(3.70) 



(3.71) 
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where the one-loop functions qn and of eq. (3.66) can be obtained directly from 
(3.48) with (3.53) and (3.54), while 



Re {SRD = 



(3.72) 



where Qn and of eq. (3.65) can be obtained directly from eq. (3.31) with (3.35) 
and (3.36). Adding (3.71) and (3.72) and using the relations Pn + hn — Qn + Qn a-nd 
Pli + K^Qn + gl we recover (3.68). 

We emphasize that qn and have formal power series expansions in terms of 
l/tj or apT(/i|). The leading term in this expansion, which is a valid approximation 
at large enough yU/ can be obtained by replacing apT(/c^) with apT(/x|) inside the 
integrals in eqs. (3.49) and (3.50), 



(3.73) 



9 \ n 



d(") 1 

BVlog-2- 



B 



(") 

R 



sine irn H (sine nn — cos Trn) 



n 



On the other hand (?„ and (?^^^ do not have similar expansions in l/tj, since they are 
proportional to a power of A'^/fij. Substituting Qn and in (3.72) one obtains a 
cutoff independent result for Re {SR^^'^}, 



Re {<5i?r}k.,oop = 



5j)logJ^ + C 



cos Tcn 



B 



in) 
R 



TT sm Tin 



(3.74) 



In general, if /// » A, 

Re{<5i?r-}|«|Re{i?-4|, (3.75) 

since Re{(5i?f^} behaves as (AVQ^)" = {n]/Q^y' (AV/xf)" ^^ile Re{i?fj,^„} be- 
haves as {ixyQ'^Y apT{fij). The only exception to (3.75) is an analytic term (integer 
n with no logarithmic terms: b]^^ — 0) where Re vanishes identically. Note, 

on the other hand, that if the leading n is half-integer and there are no logarithmic 
terms (-B^"* = 0) then Re {SR^^"^} vanishes identically. In this case i?APT coincides 
with the principal value Borel sum, up to some sub-leading power corrections which 
can be ignored if is large enough. 

3.4 Generalization to two-loop running coupling 

In the "skeleton expansion" approach (2.4) all the diagrams that correspond to dress- 
ing a single gluon are formally contained in the first term Dq^. Identifying the 
coupling in (2.5) or (3.2) with the one-loop coupling, thus amounts to some ap- 
proximation already at the level of the leading term in this expansion. In order to 
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improve this approximation, or at least to have some information about its accuracy, 
it is important to calculate the integral also with the two-loop running coupling. 

Since the "skeleton expansion" is not systematically defined, the justification of 
(2.5) or (3.2) with the particular function ^^(k'^ /Q"^) or !Fr{ii'^ /Q"^) is based on the 
large Nf limit and the "Naive Non-Abelianization" procedure. In this approximation 
the coupling is strictly one-loop, and so using eq. (2.5) or (3.2) beyond one-loop is 
not really justified. Still, having in mind the picture of the "skeleton expansion" , 
we regard the running coupling Opt inside the integral as an all-order coupling and 
eventually as a (non-perturbative) infrared regular coupling. The first stage is to 
promote it from the one-loop to the two-loop level. 

Technically, performing the infrared cutoff regularization with a two-loop run- 
ning coupling is more involved, but as we shall see in this section, it is achievable. 
One possibility is to use the exact expression [35] for the renormalization scheme 
invariant Borel transform a{z) of eq. (3.17), but the resulting integrals are not obvi- 
ously tractable. Alternatively, we use the observation that for any coupling which has 
only* a space-like Landau cut ending at k'^ — —Is? the following dispersion relation 
is satisfied 

^2^fc2 Ppt(/^')- (3-76) 

This holds in particular for the two-loop coupling with (3i > 0. It then follows [40] 
from eqs. (3.12) and (3.11) that 

(5aAPT(A;^) = / ^ -pri^ Ppt(/^^)- (3.77) 

i-A2 + k'^ 

We therefore find that the bnS in eq. (3.13) arc given by 




(3.78) 



where we exhibited the phase arising from the negative integration range. As men- 
tioned above, for integer n, the 6„'s in eq. (3.13) and thus also in eq. (3.78) coincide 
with those of eq. (3.32). We assume that the identity between (3.78) and (3.32) holds 
also for non-integer n. 

Consider now the case of a pure power contribution (i.e. no logs: 5^"^ = 0) to 
Tr. Then eq. (3.31) gives 

/ A 2 \ " 

5K'^ = -4") 6„ j (3.79) 
On the other hand, eq. (3.48) gives 

KL = -4"^ [0) " sine nn /„ (nj/^') (3.80) 
* Specifically we require absence of complex Landau singularities. 
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It was shown in [39, 42] that if Cpx satisfies the two-loop renormahzation group 
equation 

= — Po f*PT ~ Pi '^PT) (3.81) 



dlnA;2 

the standard Borel representation of is 

9i 



/. = /; exp (--i) fiiA-i^a - /; exp (3.82) 

\ ZnJ \ Z„J 



where aj = apT(/x|), 



Sn = = (3.83) 



and 

^-^ + |. (3.84) 

a-i a-i Po 

The integral in eq. (3.82) can be expressed in terms of the incomplete gamma func- 
tion, 

In = -zn (-|)'" exp (-|) r -|) (3.85) 

where 

r(z,a;) = y 7^'^"*- (3.86) 

Let us first compute the (ambiguous) imaginary part of /„. It is convenient to 
write 

\ ZnJ \ Zn) 

^ It^ + (3.87) 
where only contributes to the imaginary part. We have^ [41] 



/r^ = -^n(-^j exp (^-^j r(-5„), (3.88) 





Zn \ 








\ diJ 



and therefore 



where the identity 

r{-5n) = wtItt (3.90) 

SVa-Kdn r(l + 5n) 

'^This simple form can be obtained from eq. (3.85) by putting to zero the second argument of 
the incomplete gamma function. 
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has been used in the second hne. 

Integrating the two-loop renormahzation group equation eq. (3.81) gives 



..-.(A 



exp = C e-^" (3.91) 



where A is the tip of the Landau cut. From (3.80), (3.89) and (3.91) it follows that 

Im = ±4") sine nn n e"'" i^-^j (3.92) 

We can now compute 6„ of eq. (3.78) observing (see sec. 3.3) that this imaginary 
part must cancel with the one coming from SRa.pt, namely (eq. (3.79)) 



Im{5it:r} = -4"^Im{U(^) ■ (3.93) 



Prom (3.78), which is assumed* to be valid also for non- integer n, we obtain 

Im{U-±sin7rn|_°^^^ (^-^^ Pp,(a^^). (3.94) 
Comparing eq. (3.92) with eq. (3.93), the cancellation condition gives 



where the sign is determined by the observation that pp^di^) is negative also when 
/i^ is negative, which implies 

PO 1 (1 + On) 

Taking the real part of SRapt (eq. (3.79)) and using eq. (3.96) we get 

Remr}L^ = [q,) oosnn St e-\ (3.97) 

which generalizes eq. (3.74) with B^^^ = to two-loops. Just like in the one-loop 
case, when n is half-integer (3.97) vanishes and then i^APT coincides with the principal 
value Borel sum. 

To compute Ai? (eq. (3.69)), we next need Re{i?fj^}. We have 

(2 \ " 
sincTrnReUJ (3.98) 

*What we actually use here is only the assumption that the phase of 6„ is e^"^". 
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with (eq. (3.87)) Re{/„} = 4") + Re{/W|. But 



where 



hence 



j(z,x)= —fe-* = r(z)-r(z,x), (3.100) 
Jo t 



Re {/„} = Zn (-^) 7 (-^n, exp (^-^ 

- Zn cos 7TdnT{-6n) (^y) (3.101) 



where we also used eq. (3.88). It follows that 



+ ^^(^J cos7r5„r(-5„)5^e-^" (3.102) 

where we used eq. (3.91) to simplify the fj^j independent part. We checked that taking 
the limit Pi — in eq. (3.102), the singularities in the two terms cancel and the 
one-loop result of eq. (3.71) with (3.66) is recovered. Combining (3.102) and (3.97) 
we end up with 

Ai^nU, = -^-^ ( ^ ) (--) 7(-5n,-^jexp(^-^ 





l3o V<3V r(l + 5„) " V sm7r5, 

where the identity (3.90) has been used in the second term. The extension of these 
results to the more general case where a log term is present in the small gluon mass 
expansion of J^r can be dealt with using the d/dn trick. 

3.5 Comparing different regularizations 

It follows from eq. (3.70) that and i?pT|pv differ by infrared renormalon power 
terms contained in Re {Rfj^}, which are related to the non-analytic terms [24, 25, 46, 
8] in the small /i^ expansion of Tr. This implies that these two regularizations are in 
fact equivalent for the analysis of power corrections of infrared origin, as we shall see 
in the next section. The same statement can be made for the effective regularization 
obtained by truncating the perturbative series in a given renormalization scheme at 
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the minimum term, which turns out to be numerically close to the principal value 
regularization''" (see sec. 5.4 and fig. 7 and 8). 

On the other hand, an arbitrary regularization may differ from the principal value 
Borel sum by power terms of ultraviolet origin which arise from analytic terms in the 
small expansion of J-'r. Examples of the latter kind are the APT regularization, 
where the contribution of analytic terms to 5i?APT is apparent in (3.74) and (3.97), 
and the cutoff regularization in the Minkowskian representation R^, eq. (3.21). This 
property makes these regularizations inconvenient for the analysis of infrared power 
corrections. We stress, however, that when the leading power corrections are of the 

1/Q type (assuming n = ^ and -B^^ = 0) the APT regularization becomes quite 
convenient since it coincides (sec the end of sec. 3.3) with the principal value Borel 
sum up to some sub-leading power corrections which can be usually ignored. 



4. Power corrections 

Let us assume that the perturbative analysis of the "leading skeleton" term i?" in 
eq. (3.2) indicates a leading renormalon at z = Zn = n/^Q in the Borel plane, which 
implies that the Borel sum ambiguity is (9(1/Q^"). Since the full (non-perturbative) 
QCD result should be unambiguous, it differs at large from a generic regulariza- 
tion of the perturbative sum {R^'^) by a 0(l/(5^") power correction (see refs. [43, 44]). 
In principle, the non-perturbative corrections should be calculable in QCD. In the 
absence of such a calculation one will naturally attempt to fit experimental data by 

(4.1) 

neglecting, for simplicity, possible sub-leading power corrections. In (4.1) both the 
regularized sum of perturbation theory R^'^ and the fitted coefficient A,, depend on 
the regularization method r used. Their sum should of course be independent of r 
if the fit is successful. If two regularization methods r and r' differ by a (known) 
infrared power term Xry/Q'^"' they are equivalent, in the sense that 

Xf' — Xr -\- Xf-y (4-2) 

i.e. going from one regularization to the other amounts to a straightforward redef- 
inition of the power term coefficient. In particular we have seen that this holds for 
the principal value Borel sum i?PT|pv and the momentum cutoff regularization R^. 
Nevertheless, the latter allows for a possible physical interpretation of the power 
corrections in terms of an infrared finite coupling [7, 8]. 



'^Truncation at the minimal term is a priori a scheme dependent procedure, but we find that 
this scheme dependence is small. 
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In the approach of [7, 8] (see also [26]) one assumes the existence of a non- 
perturbative couphng 

a(A;2) = ap^(A;2) + Sa{k^) (4.3) 

regular in the infrared region. As opposed to Cpx, the non-perturbative coupling a is 
assumed to satisfy the dispersion relation* 

where /o(/U^) is the discontinuity of a, defined similarly to its perturbative part (3.3), 
and aeff(/i^) is defined by 

^p{^% (4.5) 



In the framework of the "skeleton expansion" such a non-perturbative extension 
appears quite natural: if the coupling is regular in the infrared each term Di in (2.4) 
is well defined. We shall therefore assume here that it is the "skeleton coupling" 
which plays the role of the infrared finite coupling of [7, 8]. Let us consider, as 
before, only the first term in this expansion, the equivalent of (2.5), 

D{Q') = [^a{k')^n{kVQ'). (4.6) 

One can write [16] 

+ /rl? «-(^') ^Dik-'/Q') + r^S-aik') ^nik'/Q') 

Ji?, Ju?, k'^ 



= AH(g') + DlliQ') + 6D,,{Q') (4.7) 

where is for the moment an arbitrary infrared cutoff. Following [7, 8, 26], we shall 
assume one can choose /ij in such a way that at scales above /ii the full coupling 
a is well approximated by its perturbative piece CpT. One can then neglect the last 
ultraviolet piece in eq. (4.7) 



^DM^)-rDZm- ' (4.8) 



■t While the APT coupling (3.11) has all the assumed properties for the non-perturbative coupling 
a, we do not imply that it is the correct model. 
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The infrared cutoff regularization thus appears naturally in the present frame- 
work. The power corrections arise only from the infrared piece Dir(Q^)^. This piece 
yields, for large Q^, non-perturbative "long distance" power contributions which cor- 
respond to the standard condensates for observables that admit an operator product 
expansion. If the Feynman diagram kernel ^oik'^/Q'^) is O [(/c^/Q^)"] at small A;^, 
this piece contributes an 0[(A^/(5^)"] term related to a dimension n condensate, 
with the normalization given by a small virtuality moment of the infrared regular 
coupling a{k'^) (sec cq. (4.14) below). 

The generalization of this approach to (inclusive enough) Minkowskian quantities 
has been given in [8]. One simply extends eq. (3.2) to the full non-perturbative 
coupling, to obtain at the single gluon exchange level 



Jo 



IJ, 



1^' 



The analogue of eq. (4.7) is [45, 36, 16] 



(4.9) 



(4.10) 



with 



(4.11) 



and /?u^,(Q^) defined in eq. (3.39,3.40). Assuming again that the ultraviolet piece 
dRxjy can be neglected, we end up with 



R{Q') ^ RUQ') + RlliQ' 



(4.12) 



In practice [7, 8], one expands ^R{k^/Q^) at small k'^ and obtains an approximation 
to RiYiiQ^) based on the leading power correction. Consider for example the case of 
a leading renormalon at a half integer n in (3.27) with no logarithms {b'jI^ = 0). 
Then ^R{k'^/Q'^) in eq. (3.41) is given by (see eq. (3.47)) 



(„) sm nn 

TT 



(4.13) 



and 



n ■ 



dk"^ 




d{e) 



(4.14) 



where the integral in the square brackets is a specific moment of the ^inwersa/infrared 
regular "skeleton cou]:)ling" which serves as a non-perturbative parameter and the 



§See [16] for a discussion of the more general case where the contribution of 6Duy is kept. 
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observable dependent coefficient in front comes out of the perturbative calculation 
of ^(/iVQ')- 

In this approach, the infrared cutoff fij itself acquires some physical meaning. It 
is bound to be small enough such that the approximation of Q^) by the leading 

order term in the small expansion will be valid, allowing e.g. the parametrization 
of in the form (4.14). On the other hand, /ii should be large enough such that 
a will be approximated by the perturbative coupling Opx above This would 
ensure that i?^v is under perturbative control. Moreover, the universality of the 
"skeleton coupling" a requires that for different observables a common /ij is chosen 
consistently with the requirements above. Such a choice would allow to compare the 
non-pcrturbative parameters obtained by fitting experimental data with (4.12) for 
different observables that share the same leading rcnormalon behavior. Note that it is 
possible that the data is well fitted by what appears to be an entirely "perturbative" , 
but regularized, ansatz such as i?pT|pv, or even i^^v with an unrealistically low choice 
of /ij (such as III — A!). However, this would not deter the alternative interpretation 
of the same data in terms of eq. (4.12), but this time with a more "realistic" larger 
value of Hi: we indeed saw above that the fit result cannot depend on the choice of 
Hi. On the other hand, the arbitrariness of the regularization implies one cannot fix 
the correct "physical" /ii studying of a single observable. 



5. Application: average thrust 

As an example of the method proposed we analyze here a specific observable, the 
average thrust (T) in e'^e~ annihilation. The thrust characterizes how "pencil like" 
the event is. It is defined as 

(5.1) 

Si \Pi\ 

where i runs over all the particles in the final state, pi are the 3-momenta of the 
particles and nx is the thrust axis which is defined for a given event such that T is 
maximized. For a "pencil like" 2-jet event, T approaches 1. It is therefore natural 
to define t = 1 — T, such that t vanishes in this limit. 

The definition (5.1) guarantees that T does not change due to emission of ex- 
tremely soft gluons, i.e. it is infrared safe. In addition T does not change due to 
a collinear split of a particle, i.e. it is collinear safe. These properties suggest [50] 
that the thrust distribution and the average thrust can be calculated in perturbative 
QCD from parton momenta and compared with experimental measurements where 
the thrust for a given event is obtained from hadron momenta. The gap between par- 
tons and hadrons may result in a non-perturbative modification of the perturbative 
result due to confinement effects. 
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Like other event shape parameters, measurements of the average thrust do not 
agree [3] with the next-to-leading order perturbative calculation [2, 51] 



(5.2) 



where Cp 



2Nc 



I and 



to = 1.5776 
ti = 23.7405 - 1.689 



(5.3) 



The significant discrepancy is shown in fig. 1. We note that the experimental data 
points (there are 44 data points altogether) are fairly consistent between different 
experiments, and show a stronger variation with Q. 
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Figure 1: Average of 1 — thrust as a function of the center of mass energy Q, according to 
the available experimental data [52] and the leading order (LO) and next-to-leading order 
(NLO) perturbative QCD calculation in the MS scheme with fiji = Q, given (Mz) = 
0.117. 



There are theoretical indications that perturbative corrections at the next-to- 
next-to-leading order and beyond are large: the next-to- leading order term is quite 
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a significant correction with respect to the leading order and there is a large renor- 
malization scale dependence. However, the usual explanation of the discrepancy 
between the data and the next-to-leading order calculation is that hadronization re- 
lated power corrections are important. It is intuitively clear from the definition of 
the thrust that emission of soft gluons implies a change in the thrust which is linear 
in the gluon momentum, and thus power like effects that fall asl/Q can appear. Fits 
to the experimental data [3] confirm the existence oi 1/Q power corrections. Event 
shape variables are quite special in having such strong power corrections: other 
QCD observables, like the ones for which higher twist terms can be analyzed using 
an operator product expansion, usually have power corrections that fall as or 
faster. 

Traditionally, extraction of ag from experimental data is based on Monte-Carlo 
simulations that effectively generate "hadronization" power corrections, which are 
added to the next-to- leading order perturbative expression (5.2). In the last 5 years 
there have been several theoretical attempts to understand better the subject of 
power corrections in event shape variables. These includes analysis of the interplay 
between renormalons and power corrections [4, 5, 6, 9] and renormalon inspired ap- 
proaches to parametrize power corrections, based on either a universal infrared regu- 
lar physical coupling [7, 8] or a less restrictive observable-dependent non-perturbative 
shape function [10]. 

Since the average thrust is a priori expected to have both large perturbative 
corrections and significant power-like corrections, we find this observable quite ap- 
propriate to serve as an example of our approach. 

5.1 The characteristic function for the average thrust 

The observable dependent ingredient in the renormalon resummation program is the 
characteristic function r{^^ / Q"^) ■, which is defined by the leading order perturbative 
result for a gluon of mass /i.^. In this section we shall calculate the characteristic 
function for the average thrust, Tt^I^^ IQ'^)- 

As mentioned in sec. 3, Minkowskian quantities usually have two different ana- 
lytic functions: ^^(-^(e) for e = n^jQ^ < 1 and J^^+\e) for e > 1. The reason is that 
the former includes both real and virtual gluon diagrams while the latter includes 
only virtual gluon diagrams. In case of event shape variables like the thrust, virtual 
corrections do not contribute at one-loop, and so J-'^\e) = 0. What remains to 
calculate is JF^, -"(e) which is entirely due to real gluon emission. Let us simplify 
the notation and define J-'{e) = J-'^ \e). The "gluon mass" integral (3.2) is now 
performed up to /^^ = 



which implies that there are strictly no ultraviolet renormalons in this case: the large 
order behavior of Rapt{Q^) when expanded in some scheme is determined just by the 




(5.4) 
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leading non-analytic terms in the small e expansion of J?^(e), which are the leading 
infrared renormalons. The absence of ultraviolet renormalons is a direct consequence 
of the absence of virtual corrections at the order considered. Higher terms in the 
"skeleton expansion" may give rise to ultraviolet renormalons in the full perturbative 
series of (t). 

The characteristic function T{e) was calculated numerically in [8] , and its leading 
term in the small e expansion was obtained there analytically, j^{e) — 4-^, which 
indeed indicates 1/Q corrections. We shall repeat the calculation here. We find it 
important to have, if not an analytic expression for J-'{e), then at least several leading 
terms in its asymptotic expansion for small e. The sub-leading terms are required 
to verify the convergence of the power correction series (we shall see that in practice 
only 1/Q terms are important). 

In order to explain how J-'{e) is computed we first briefiy review the kinematics 
and the calculation of the thrust for three partons in the final state. Let us denote the 
primary photon 4-momentum by Q, the 4-momenta of the quark and anti-quark by 
Pi andp2 (Pi = 0), and the 4-momenta of the "massive gluon" by p3 (pf — ii^ = eQ"^). 
Following [2] we define yij = (pi+Pj)'^ / Q'^ , which implies i/ij > 0, and Xi = 2Q ■ Pi/Q^. 

It follows from energy-momentum conservation that 

Xi + X2 + Xs^2 (5.5) 

and from the assumed virtualities of the particles that 

yi2 + X3 - e = 1 
yi3 + X2 = I 

1/23 +Xi = 1. (5.6) 

In the center of mass frame, where Q — (Q, 0, 0, 0), Xi becomes the energy fraction 
of the i parton, Xi — It then follows that 




i = 1,2 



\ps\ = ^^j4-^e. (5.7) 
The phase space limitations are the following: 

a) Hard gluon limit, 

Xi + X2>l-e, (5.8) 
which follows from the condition yi2 > together with (5.5) and (5.6). 

b) Soft gluon limit, 

(1 - xi)(l - X2) > e, (5.9) 
which is obtained in the center of mass frame, using (5.7), from the condition 

IP2I = |pi +P3I < \pi\ + IpsI- 
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The phase space hmitations are shown in fig. 2 for the case e = 0.1. The region 
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Figure 2: Phase-space for the emission of a virtual gluon with = eQ^ = Q.IQ^ in 
the plane of the quark and anti-quark energy fractions (xi.2). Continuous lines represent 
phase-space limits: the upper (curved) line corresponds to the softest gluons (5.9) while the 
lower (linear) line corresponds to the hardest (5.8). Dashed lines represent the separation of 
phase-space according to which particle carries the largest momentum and thus determines 
the thrust axis (cf. eq. (5.13)): in the upper left region T = X2, in the upper right region 
T = xi and in the lower region T = y^Xg — 4e. 

near the outer (curved) line corresponding to (5.9) represents soft gluons: for a given 
"gluon mass" eQ^ and quark energy fraction (xi), the smallest possible gluon energy 
is obtained when the inequality (5.9) is saturated. The region near the inner (linear) 
line corresponding to (5.8) represents hard gluons with maximal energy X3 = 1 + e. 
Note that as e increases the relevant phase space shrinks and approaches the region 
of small Xi and X2 (the lower left corner in fig. 2). 

The characteristic function J^(e) is obtained from the following integral over 
phase space [8], 




phase space 




(5.10) 
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where ^ a is the squared tree level matrix element for the production of a quark, 
an anti-quark and a gluon of "mass" fi^ = eQ^ and 

M(xux^,.)= (1 _,^) -(13^ -(13^- (5-11) 

The last ingredient for the calculation of T{e) is the expression for the thrust. 
For two particles in the final state, the thrust axis ut coincides with the line along 
which they move and T — 1. For three particles in the final state, it coincides with 
the direction of the particle carrying the largest momentum, \pi\. By momentum 
conservation in the center of mass frame, the two other particles have the sum of 
momenta Pj+Pk = —Pi- The numerator in (5.1) then equals 2\pi\. For three massless 
particles the denominator in (5.1) equals \pi\ + \p2\ + jpsl = Ei + E2 + = Q. 
Thus, using (5.7), we have T = Xi. For a massive gluon, however, the denominator is 



^{xi+X2 + y xl — 4e) 7^ Q. On the other hand, since the "massive gluon" dissociates 
into massless partons one should actually calculate the thrust taking into account the 
final partons. In this case the denominator is Q. Note that the numerator remains 
the same whether or not one takes into account the gluon dissociation, provided that 
all the partons produced end up in the same hemisphere as the parent gluon - an 
assumption to which we return below. In conclusion we find that, due to the change 
in the denominator of (5.1), by giving mass to the gluon (in order to represent its 
dissociation through the dispersive approach) we unwillingly change the calculated 
value of the thrust. To solve this difficulty it has been suggested [8] to modify the 
definition of the thrust such that the normalization will be with respect to the sum 
of energies (Q) , 

rp ^ E» \Pi ■ ftrl ^ Ei \Pi ■ firl , . 

EiEi Q - ^ - ' 

There is no difference between (5.12) and (5.1) so long as only massless particles 
are produced. However, for the theoretical computation with a "massive gluon" it 
is important to use (5.12) which guarantees that the same value of the thrust is 
obtained with a "massive gluon" as with the massless products of its dissociation^. 
The final result for the thrust, using (5.12), is 

t{xi,X2, e) = min |l - Xi , 1 - X2 , 1 - \J {2 - xi - ^2)^ - 4e| . (5.13) 

Fig. 2 shows the separation of phase space to regions where each of the three particles 

has the largest momentum and thus determines the thrust axis. 

Let us now return to the more delicate problem, namely the assumption we made 
concerning the decay products of the gluon. This assumption is absolutely necessary 



^Thc original definition of the thrust, eq. (5.1), with a "massive gluon" does not comply with 
this requirement. It would lead to a different result for ^(e) and thus to a different normalization 
of the power corrections (see [47, 8]). 
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in order to keep the same value of the thrust when referring to the "massive gluon" 
itself as when referring to the massless products of its dissociation. In the first 
case the relevant term in the numerator of (5.12) is |p3 ■ utI while in the second it 
is larger: |p2 ■ utI + \pr • ^t|, where ps = Pl + Pr and {pr) corresponds to the 
sum of momenta of the particles that originate in the gluon and end up in the left 
(right) hemisphere. It is a priori not clear whether this assumption is justified: 
kinematic considerations alone do not exclude the possibility of dissociation into 
opposite hemispheres and in fact when the gluon is close to the transverse direction, 
it is quite plausible that it would dissociate this way. The "massive gluon" approach 
is justified [24, 25, 8] for completely inclusive quantities: then dressing the gluon 
or taking into account its dissociation amounts exactly to building up the running 
coupling. We see that the thrust is not inclusive enough. We can still ask, however, 
how large is the error we introduce using the inclusive "massive gluon" approach 
instead of taking into account separately the contribution of the decay products of 
the gluon. The problem of non-inclusiveness was first raised in [5] in the framework 
of renormalon resummation in the large Nf limit where the terms which prohibit an 
inclusive treatment were identified and evaluated. We shall return to this issue in 
the next section. 

Let us now proceed with the calculation of J^(e) for the thrust, based on the 
definition (5.10), the squared matrix element (5.11) the expression for the thrust 
(5.13) and the phase space limitations (5.8) and (5.9) which determine the integration 
range. In order to evaluate the integral (5.10) one has to treat separately each of the 
three regions of phase space (see fig. 2): the integrand in each of them is different, 
as imphed by (5.13). To perform the integrals it is useful to change the integration 
variables to the following: zi — xi + X2 + e — 1 and Z2 — xi — X2, which fit better 
the phase space limitations. Most eventual integrals can be performed analytically 
but the resulting functions are complicated. Instead, we calculated the integral 
numerically for any < e < 1 and in addition obtained an asymptotic expansion of 
J-{e) for small e. 




-—+77^ + 8 ln(3) ln(2) - ^ ln(3) + 4 dilog(4) + 4 dilog(3) 
18 4 

1 r / M 160 3 

8e2 + [4+ 121n(3)] e - — 



(5.14) 



33 



5 

;e2 + 



-0.3061 + 5.3333 In 



-llln(^) + 3.4709 



The agreement between the asymptotic expansion (5.14) and the numerical calcula- 
tion is shown in fig. 3. 
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Figure 3: The characteristic function J^(e) for the average thrust (5.10) as a function of 
logio(e)) where //^ = eQ^ is the "gluon mass". J^(e) is represented by a thick black hne 
reaching zero in the hmit e — > 1. It is compared with its asymptotic expansion at smah e 
(5.14), for n = ^, 1, |, 2, |, where n is the highest order term ©(e") taken into account in 
each approximating curve. 

As we saw in the previous sections, the leading terms in the small e expansion of 
jF(e) are required for the calculation of the difference between different regularizations 
of the perturbative sum and eventually for the parametrization of power corrections. 
Each term in eq. (5.14) has the form of eq. (3.27) and so the numerical coefficients 
B^^^ and 

can be immediately identified. The first non-analytic term in (5.14) is a 
square root one, corresponding to a 1/Q infrared power correction. The coefficient of 
this term agrees with previous calculations [8]: ^(e) = 4:y/e. As explained there, this 
term arises from the limit of phase space (5.9) corresponding to the softest gluons. 
A new observation is that there are no l/Q^ infrared power corrections in this 

3 

"leading skeleton" approximation. The next non-analytic term in (5.14) is asso- 
ciated with infrared power corrections. As a result the apparent convergence 
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of the discontinuity function ^{k'^/Q'^) 



TT 




8 - 17.7778 



(5.15) 



is better than that of J^{e). In (5.15) the second term becomes about 20% of the first 
around fc^ ~ 0.1 Q^- Another observation is that the asymptotic expansion (5.14) 
does not contain double logarithmic terms up to the order considered. 

Finally, taking the logarithmic derivative we obtain J-'{e) = —edJ-'/de, which 
agrees with the numerical results of [8]. In fig. 4 one can see that there is just a small 
contribution to the characteristic function, and thus also to the value of the average 
thrust, from the kinematic configurations where the gluon momentum is the largest. 
Moreover, this contribution to J^(e) has some significance only for fairly large "gluon 
mass" e ^ 10~^, and so its effect on the large order behavior of the perturbative 
series and the related power corrections is negligible. 



ID 

I 




Figure 4: The derivative of the characteristic function for the average thrust J-{e) as 
a function of Iog]^o(^)) where /i^ = eQ^ is the "gluon mass". The separate contributions 
from kinematic configurations where one of the quarks or the gluon carries the largest 
momentum are shown as well. The dashed line represents the leading (n = ^) term in the 
small e expansion of this function, .F(e) ~ 4A/e. 
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5.2 Does the renormalon resummation program apply to non-inclusive 
quantities? 

The "skeleton expansion" approach and thus our renormalon resummation program 
are based on the assumption that all the higher order diagrams related to the running 
of the coupling contribute inclusively to the observable. This assumption does not 
hold in the case of event shape variables like thrust. In the previous section we 
calculated the thrust characteristic function using a "massive gluon", which should 
represent after integration with respect to the discontinuity of the running coupling, 
the partons to which the gluon dissociates. We encountered and eventually ignored 
the non-inclusive nature of the thrust: as explained there, if the decay products 
of the gluon end up in different hemispheres the average thrust value obtained in 
the "massive gluon" calculation is different (and in fact always lower) than it is in 
reality. Still, we intuitively expect that most decays are roughly collinear and so 
for the major part of the phase space decay into opposite hemispheres is not very 
likely. We therefore conjecture that the inclusive resummation approach can provide 
a reasonable approximation to the full higher order corrections. In the following 
we give a quantitative argument at the next-to-leading order level in favor of this 
conjecture. 

In order to estimate the error we are making by the inclusive treatment, let us 
compare the two type of expansions we have: the ordinary perturbative expansion, 
e.g. eq. (5.2) in the MS scheme, and the conjectured "skeleton expansion", i.e. the 
equivalent of (2.4) and (3.1) for the average thrust 

(^)PT iQ^) = ^ [[sr{Q^) + Sr{Q^) + ...)+ Non - Skeleton] (5.16) 

where similarly to R\'^ in eq. (3.1) S^'^ is the leading "dressed skeleton" term, nor- 
mahzed as -S'o'^ = a-\-0{a'^), S'f'^ is the second "dressed skeleton" term which starts 
at order a? and so on. The Non — Skeleton piece allows for additional perturbative 

terms of order or higher which do not fit into a "skeleton expansion" due to the 
non-inclusive nature of the observable. The "leading skeleton" term S'g'^ is given by 
i?APT of eq. (3.2), up to a regularization related power term which we now ignore. 
Since both S'f^ and the Non — Skeleton terms start at order a^, we can approximate 

(^)p. m by 

with 

5nlo = W + • • • (5.18) 

where in general the coefficient ii contains contributions from both and the 
Non — Skeleton terms. To compare the two expansions, let us expand the time-like 
couphng within the first "skeleton" term -Rapt (3.2) in terms of Cms ■ Note that at 
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a difference with sec. 2, it is the Minkowskian representation that we start with. 
Taking the one-loop form (3.5) of the coupling we obtain from (3.6) the following 
expansion in (Q^) 

,2 \ "1 



Klif^') = «Ms {Q') + \d^-{c^ + \og^]f3o}a|^, (Q^) + 



di^ - 2di fci+log^j/3o + 



(5.19) 



r 



where the first two terms are actually valid beyond the one-loop approximation of 
the coupling and also coincide with the corresponding terms in the expansion of the 
space-like coupling (2.7). Inserting (5.19) into (3.2) yields 



with 



Sr{Q')^toajrs{Q') + i?a|s(Q') +^2«ls(Q') + 



*o — /o 

i? = c^i/o-(ci/o + /i)/5o 
4 = di^ fo - 2di{cifo + f^)(3o+ [ci 



(5.20) 
(5.21) 



77 tt' /o + 2ci/i + /2 



and so on. Note that at the next-to-next-to-leading order (ig) we recover the char- 
acteristic TT^ terms which appear in perturbative expansions of Minkowskian observ- 
ables. The constants /j are the log-moments of the characteristic function (compare 
with (2.10)) 

dfi^ 



fi 



log^) ^(/.VQ^). 



(5.22) 



Using the numerical result for the characteristic function of the average thrust one 
can obtain fi to any arbitrary order. The first values are 



/o 
/i 

/2 
/3 
/4 

h 
h 



1.577602558 

-7.176762311 

42.11235577 

-307.991760 

2736.14010 

-28923.5429 

357358.9993 



(5.23) 



With these coefficients at hand one can construct a power series approximation!' 
(5.20) to the first "skeleton" term S^'^ of eq. (3.2) to any arbitrary order, provided 

"This expansion is an asymptotic one, badly affected by infrared renormalons. Note that tfie 
explicit sign oscillation in /j cancels against the sign oscillation in eq. (5.21). Note also that the 
fast growth, that eventually becomes factorial, is already apparent in (5.23). This expansion will 
be discussed in sec. 5.4. 
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one specifies the "skeleton coupling" OpT, namely the parameters di and Ci. As 
discussed at the end of sec. 2, we use in this work several different schemes for the 
"skeleton coupling" ap-p. In the Abelian limit, apT should coincide with the V-scheme 
coupling and so Ci = — |. We can thus determine the (3q dependent piece in i° in 
eq. (5.21) 

-(ci/o + /i)/3o = 9.8061 /3o. (5.24) 

We now have all the ingredients for the comparison up to next-to-leading order 
between the "skeleton expansion" (5.17) and the standard expansion (5.2). For the 
first we use (5.20) and obtain 

(^)PT = ^ {^0 ajTs iQ') + {ii + t?) airs iQ') + O (a|, (Q^)) } (5.25) 



2 
Cf 



{fo Oms {Q') + [ii + d,fo- (ci /o + /i) Po] {Q') + O (alg {Q')) } 



By construction the leading order is the same. The comparison at the next-to-leading 
order gives 

ii + difo + 9.8061 Po = -4.128 + 10.134 Pq (5.26) 

where the l.h.s. corresponds to the "skeleton expansion" coefficients ti + of 
eq. (5.25) and the r.h.s. to the standard expansion coefficient ti of eq. (5.3), with 
the Nf dependence expressed in terms of /3o (2.3). 

For an inclusive quantity, where we assume that the "skeleton expansion" exists 
(i.e. the Non — Skeleton terms in (5.16) are absent), the Nf dependence at the next- 
to-leading order comes only from diagrams that arc related to the running of the 
coupling. In such a case the entire Pq dependent piece in the next-to-leading order 
coefficient is accounted for by the leading term 5*0^ in the "skeleton expansion" , and 
the remaining coefficient ti in (5.17), which coincides now with the normalization of 
the sub-leading "skeleton" S^'^ = iio^ + • ■ ■, should be free of Pq. For the thrust, 
which is non-inclusive with respect to the decay products of the gluon, this does not 
hold. However, we find that the difference between the term linear in Po in (5.20) 
9.8061 Po (l.h.s. in (5.26)) and in full next-to-leading QCD coefficient 10.134 Po (r.h.s. 
in (5.26)) is quite small: it is about 3 percent. This finding gives place to hope that 
the Non — Skeleton terms in (5.16) are small and thus the inclusive treatment is after 
all a good approximation for the resummation of a certain class of diagrams. 

The observation that for a non-inclusive quantity the Nf dependence of the next- 
to-leading order coefficient cannot be explained in terms of the running coupling 
also implies that the usual motivation for BLM scale fixing and the Naive Non- 
Abehanization procedure** does not hold. As opposed to the inclusive case discussed 
in sec. 2, the BLM scale computed from the full next-to- leading order coefficient does 



**The ambiguity of the Naive Non-Abelianization procedure for non-inclusive quantities was 
pointed out in [48]. 
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not coincide with the one of the "leading skeleton". As a result, the identification 
made following eq. (2.13), between the remaining next-to-leading order coefficient^^ 
and the BLM coefficient in the "skeleton scheme" fails: the former still contains 
some Nf dependence while the latter does not. Like our resummation program, the 
BLM procedure becomes relevant once the observation is made that the /?o dependent 
terms in the two sides of eq. (5.26) are numerically very close. 

Finally, using (5.26) we evaluate ii. If the coupling a^^ in the "leading skeleton" 
term (3.2) is the time-like coupling associated (3.4) with the "gluon bremsstrahlung" 
couphng [30] (see the end of sec. 2), where d\"'^ — 1 — then 

^-brem ^ _i_8i3 + 0.328 Po. (5.27) 
For the pinch technique coupling [22] d^'""'' = 1, and so 

^pinch ^ _5_706 + 0.328 po- (5.28) 

For the V-scheme couphng [15], dX — —2, and then 

iX = -0.973 + 0.328 Pq. (5.29) 

For Nf = 5 the coefficients are i^"''^ = -1.184, tp*" = -5.077 and = -0.344. 
These coefficients can be compared with the standard next-to-leading coefficient 
in MS (5.3) which equals ti ~ 15.296. We conclude that at least the apparent 
convergence of the suggested expansion is better than that of the standard one. The 
coefficient is extremely small in the "gluon bremsstrahlung" and V schemes. We 
shall thus choose for our phenomenological analysis two couplings: 

a) the "gluon bremsstrahlung" scheme which was used before in the analysis of 
power corrections to the thrust [7, 8] . This coupling will be quite convenient in 
practice: having a small next-to- leading coefficient (5.27), the full result should 
be close to the first term in the "skeleton expansion" . 

b) the pinch technique coupling which may be the correct physical "skeleton cou- 

pling" a. Using this coupling, with its relatively large next-to-leading coefficient 
(5.28), in addition to the "gluon bremsstrahlung" couphng, will be useful to 
measure the sensitivity of our procedure to the value of di. 

5.3 The perturbative sum vs. experimental data 

The most convenient way to calculate the perturbative sum is to use the APT formula 
(3.2). As explained in sec. 3, other regularizations of the perturbative sum, such as 

the ]:)rinci]:)al vahie Borel sum can then be obtained from Rapt- 

There ri — r° was also identified with the normalization of the "sub- leading skeleton" . 
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The ingredients required for the calculation of -Rapt are the numerical function 
J^ili'^/Q^) - T{0) we obtained in sec. 5.1 and the discontinuity of the perturbative 
coupling on the time-like axis, Ppt(a*^)- For the latter, we shall use here the one and 
two loop couplings. In the one-loop case (3.5), the expression obtained from (3.3) is 
simply 



/?0 7r2 + log 2(g)' 



In the two-loop case we use the Lambert W function representation of the coupling 
[37, 38] 



i2-loop(^ ) - 




A 1 + W^i{z) 

(5.31) 



where W{z) is the Lambert W function defined by W{z) exp [^^(-z)] = z and the 
particular branch W-i{z) is implied by asymptotic freedom: in the ultraviolet z — > 
0~ and W_i{z) — > —oo [38]. Note that in (5.31) the explicit Landau pole and the 
tip of the branch cut coincide {W-i{—l/e) = —1) and so there is only one singularity 
in the complex momentum plane, at k'^ — A^, with a cut — oo < k'^ < A^. Using the 
computer algebra program Maple, W-i{z) is readily available at any given complex 
z. It is then straightforward to obtain the time- like discontinuity P2-ioop(a*^) (3-3) 
corresponding to a2-ioop(^^)- 

As a first trial, let us calculate -Rapt based on the world average value of as, 
Q,Ms (^M^) = 0.117. We choose in the "gluon brcmsstrahlung" scheme. Taking^ 
Nf = 5 we find A^;';^" = 0.130 GcV and A^;;:^" = 0.361 GeV. We evaluate the 
perturbative sum .Rapt(Q^) in (3.2) by a numerical integration of J^{jj,'^/Q^) — J^{0) 
times either -pi.io.^W A-ioo^if^"^) ■ 

Next, we use .Rapt to calculate the principal value Borel sum, according to 
-RpT|pv = -Rapt — Re{5i?APT}, where Re{5i?APT} is evaluated for a generic term in 
the small e expansion of ^(e) (5.14) by eq. (3.74) and (3.97) in the one and two loop 
cases, respectively. At order n in the expansion, the contribution to Re{(5i?APT} is 
C)(l/Q2»). Since Re in eq. (3.74) and (3.97) vanishes identically, the leading 

contribution is at order n = 1. We note that the n = 1 term in T{e) is analytic, and 
so this contribution is not of infrared origin. We obtain 



A^ 1 A^ 

Re{5i?APT}L = Ci^ V = 8.96 

Rc {5/?apt}|., = 4'^^ r(l + 6,) e-'^ = 3.18 (5.32) 



^This point will be discussed in sec. 5.6. 
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where the numerical vahies were obtained using; = 13-^/^1 and = 17.1833. We 
find that Rc{5i?APT} is absolutely negligible, already at the lowest relevant experi- 
mental value Q = 12 GeV, where it is less than one percent of i^APT (see table 5 in 
sec. 5.5). 

The results for the average thrust (t)p^ (Q^) — ^-Rapt — ^-Rpt|pv are shown 
in fig. 5 (within the resolution of the figure i?APT and i?PT|pv could hardly be dis- 
tinguished). The first observation is that the difference between the one-loop and 
two-loop resummation results is quite small. This stability, which shall be discussed 
further in sec. 5.6, is reassuring since in our approach a^^ should actually be an all- 
order running coupling, and so its replacement by the one-loop coupling at all scales 
is not obviously justified, as already noted in sec. 3.3. 

In order to compare our results with experimental data, we should include the 
estimated contribution of the terms not included in the "leading skeleton" according 
to eqs. (5.17) and (5.18). For a concrete estimate of 5nlo we replace the arbitrary 
scheme coupling a by* the natural effective charge at hand, namely the value of the 
"leading skeleton" with the appropriate normalization 

(^NLo = iia^ ^ ii (5.33) 



and thus 



Rapt + tj 



fo 

Ra 



(5.34) 



fo J 

This expression exhausts our knowledge concerning the perturbative contribution 
to the average thrust, as it includes, in addition to the resummation of the first 
"skeleton" , the full next-to-leading order coefficient. As seen in the figure, the fine 
representing (5.34) does not deviate much from the "leading skeleton" results. This 
is due, of course, to the small ti coefficient (5.27). 

The next crucial observation in fig. 5 is that the resummed results turn out to 
be quite close to the experimental data, significantly closer than the next-to-leading 
order result in MS with fiR = Q given in eq. (5.2). As implied by the 1/Q nature of 
the leading renormalon term in the expansion (5.14) we introduce a non-perturbative 
parameter A and add an explicit power correction of the form X/Qto the perturbative 
prediction (5.34) 

(i) = ^[i?APT + 5NLo] + ^. (5.35) 

Being unable to compute A from the theory, we determine it by performing a fit of 
(5.35) to the data. The results of such a fit are summarized in table 1, where -Rapt is 
calculated with the "skeleton coupling" (assumed to be the "gluon bremsstrahlung" 
coupling) at one or two loops. 

* An alternative choice (that would be numerically close) is to replace a by the "skeleton coupling" 
ttpT at the BLM scale (5.42). 
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Figure 5: Average of 1 — thrust as a function of Q: experimental data is compared 
with naive perturbative QCD results (LO and NLO in MS with /xr = Q) and with the 
resummed perturbative series in the APT or principal value Borel sum regularizations (the 
two coincide). All the theoretical calculations are based on (Mz) = 0.117. Both one- 
loop (upper dashed) and two-loop (dot-dash) resummation results (Rapt) are presented, 
where the running coupling ap^ is in the "gluon bremsstrahlung" scheme. In the two-loop 
case we also show eq. (5.34) (lower continuous line; just below the upper dashed line) 
which includes the full 0{a1) term and eq. (5.35) (upper continuous line) which includes 
in addition a fitted X/Q term. The lower dashed line is the absolute value of the imaginary 
part of the one-loop Borel sum, which reflects the magnitude of the renormalon ambiguity. 



— PT 

a „ 

en 


A (GeV) 


/point 


one — loop 


0.36 


2.43 


two — loop 


0.16 


3.46 



Table 1: Power term (A) fit results with af^ (M^) = 0.117 based on (5.35). 
The fit results in the two-loop case^ are presented together with the perturbative 



tXhe fit results in the one-loop case (not shown in the plot) are very close to those of the two- 
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sum in fig. 5: they turn out to be quite close. We conclude that a major part 
of the discrepancy between the next-to-leading order result and the data is due 
to neglecting higher order perturbative corrections that can be resummed using the 
suggested program. This finding will be discussed in more detail in the next sections. 

A closer look at fig. 5 reveals that our theoretical results undershoot the low Q 
data points, while they overshoots at least part of the high Q data points. Indeed 
performing a two parameter fit where both the coupling and the power term are free 
can lead to better agreement with the data. This is shown in table 2. 



^PT 


«r (Mz) 


A (GeV) 


xVpoint 


one — loop 


0.111 


0.73 


1.33 


two — loop 


0.110 


0.(32 


1.35 



Table 2: Fit results for ag and A based on (5.35). 

Fig. 6 shows the perturbative summation results for a^^ (Mz) = 0.110 together 
with the best fit line (5.35) in the two-loop case. While in fig. 6 the perturbative 
sum (5.34) by itself is not as close to the data, it is clear that the main conclusion 
we drew from fig. 5 concerning the significance of the resummation holds. 

The curvature of as a function of ag and A around the minimum reflects the 
spread of the experimental results. The two parameter fit of (5.35), calculated at 
two-loop, to the entire set of data points yields an experimental error of 

a^(Mz) = 0.110 ±0.0017, (5.36) 

and 

A = 0.62 ±0.12 GeV (5.37) 

for a confidence level of 95%. 

Further statistical analysis shows that the various experiments are quite consis- 
tent and that there is no significant difference between small Q and large Q data 
points as far as our fit is concerned. For instance, if we exclude the lowest data 
points Q < 22 GeV, which seem quite spread, we find an improvement in the fit with 
a minimal x^/point = 1.15 but the corresponding values of Ug and A change just a 
little, namely: af^ (Mz) = 0.111 and A = 0.54 GeV*. If we exclude the highest data 
points Q > 172 GeV we find x^/point = 1.44 with the same central values as in (5.36) 

loop case. The difference between the best fits in the two cases is of some significance only for 
Q ^ 30 GeV and it reaches 6% at the lowest data point Q = 12 GeV. The difference at low Q 
explains the variation in in tabic 1. We further comment on the comparison between the one 
and two loop resummation results in sec. 5.6. 

*The latter central values are obtained also if we exclude the 4 data points of the Mark J 
experiment which are higher than the rest. Indeed here the fit is better: x^/point = 1.01. 
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Figure 6: Average of 1 — thrust: experimental data vs. naive and resummed QCD 
predictions, given af^ (Mz) = 0.110 - the value which yields the best fit in eq. (5.35). 
The resummation is performed with Opx in the "gluon bremsstrahlung" scheme. See the 
caption of fig. 5 for further details. 

and (5.37). The most striking evidence that there is no systematic trend in the data 
which is missed by our fit is that even if we exclude all the data points above or below 
Mz = 91.2 GeV, the best fit values are hardly affected: in the former case, having 
29 data points, we get x^/point = 1.61 with the same central values as in (5.36) 
and (5.37) and in the latter, having 20 data points, we get x^/point = 1.17 with 
af^ (Mz) = 0.111 and A = 0.53 GeV. Note, however, that the effective experimental 
error changes significantly, e.g. in the latter case, the error in a two parameter fit 
with 95% confidence level on the extracted value of as becomes ±0.011. 

5.4 Truncation of the perturbative expansion 

Let us consider now the expansion of the renormalon integral (3.2) in some renor- 
malization scheme, e.g. the expansion (5.20) in aj^{Q'^). As explained in sec. 5.2 
this series diverges at large orders due to the factorial increase of the coefficients 
induced by (infrared) renormalons. A standard procedure dealing with asymptotic 
expansions is to sum the series up to the minimal term. This can be regarded as an 
effective regularization of the all order sum. 
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As an example we analyze here the expansion (5.20) in some detail. Fig. 7 
shows the increasing order partial sums at a given center of mass energy Q = Mz 
while fig. 8 shows the results obtained when truncating the series at the minimal 
term as a function of Q. In the latter, the relevant curve is made of four distinct 
pieces according to the number of terms included in the sum: the minimal term is 
reached between the sixth and the ninth term depending on Q. Both figures show 
that truncation at the minimal term is quite close to the principal value Borel sum 
regularization. 
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Figure 7: Partial sums corresponding to the expansion of the one-loop renormalon integral 
(3.2), with cipT in the "gluon bremsstrahlung" scheme, in terms of Oms (A'/j) ^ a function 
of the truncation order, for Q = Mz = 91.2 GeV with a^(Mz) = 0.117. The lower line 
corresponds to /i^ = Q^, eq. (5.20), and the upper line to fij^ = ^bl^' ^1- (5.40). The 
square symbol is the minimal term in each expansion. The horizontal band represents 
Borel summation, where the middle line is the principal value regularization and the two 
dashed lines show the estimated renormalon ambiguity based on the imaginary part of the 
Borel sum. The symbols on the left show the experimental data points. Note that the 
small next-to-leading order correction of (5.18) is not included in the theoretical results 
shown here. 

It is clear from fig. 7 that the contribution of the 0(a|jg) and C(a^) terms is 
significant, as one could already guess based on the largeness of the 0{a^) term in 
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Figure 8: Comparison between different regularizations of the perturbative sum with 
in the "gluon bremsstrahlung" scheme for (Mz) = 0.110 as a function of the 
center of mass energy Q. The regularizations shown are: one-loop (dashed) and two- 
loop (dot-dash) principal value Borel sum (or /Japt)) truncation of the MS expansion 
(5.20) at the minimal term and the two-loop ii^y with the cutoff scales = 2 GeV and 
fij = A^^°°^ = 0.361 GeV. Note that the small next-to-leading order correction of (5.18) is 
not included in the perturbative sum in this plot. 

(5.2). The minimal term, which turns out to be close to the principal value Borel 
sum, is obtained for this center of mass energy at order 0{aj^). In (5.20), a new 
log-moment fi enters at each order in the perturbative expansion, such that the term 
(!?(a|jg ) depends on all fj with j < i. From the definition of the log-moments (5.22) 
it follows that the larger the order i, the more sensitive fi is to the small e behavior of 
jF(e). Eventually, when the asymptotic regime is reached, the added terms reflect just 
the leading term in the small fi^ expansion of the characteristic function, namely the 
leading renormalon, which is related in our case to a 1/Q power correction. Since 
the value of the added non-perturbative term X/Q is determined by a flt, it may 
not be important at which order exactly the series is truncated, provided it is close 
enough to the asymptotic regime. On the other hand, a qualitative difference exists 
between truncation of the series at the minimal term and truncation much before 
the asymptotic regime, say at the next-to-leading order. The latter would not differ 
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from a generic regularization of the pcrturbative sum just by power terms. This also 
imphes that the value of Us obtained by fitting the data with a next-to-leading order 
series plus a power term would be different than in the current approach. 

Next consider, as a pedagogical exercise, a fit to experimental data based on the 
truncated expansion (5.20) in the MS scheme, namely 



{t) 



Cp_ 
2 

CV 
2 



qpt I I r 

•^0 \o(a*') 



+ 



A 



pert 



Q 



(5.38) 



to Oms (Q') + tx airs (Q') + E til ^irs (Q') 



+ 



A 



pert 



Q 



where k is the order of truncation. The fit results are hsted in table 3 for A; = 2 
through 6. For these values of k the series is still convergent: the diverging part of 
the expansion is not reached for any Q, since the minimal term is between the sixth 
and the ninth term. The coupling a^g (Q^) in (5.38) is assumed to obey the two-loop 
renormalization group equation with Nf = 5. 



k 


«r (Mz) 


Apert (GeV) 


2 


0.128 


0.72 


3 


0.118 


0.65 


4 


0.115 


0.58 


5 


0.114 


0.50 


6 


0.114 


0.40 



Table 3: Fit results based on (5.38) using the expansion (5.20) up to order k. 

Note that nothing can be learned from the quality of the fit: it is roughly the 
same in all cases, x^/point ~ 1.3, and it is also very close to the resummation based 
fit of the previous section, x^/point ~ 1.33. The corresponding experimental error in 
the extracted value of as in table 3 based on a two parameter fit with 95% confidence 
level is ±0.0024. 

Care should be taken comparing the results for A; = 2 in table 3 with the fit 
in [7] as well as with recent experimental fits [3]. In the latter, the finite infrared 
coupling formula of refs. [7, 8, 45] is used, namely the coefficient in front of the 1/Q 
term is modified, A — > A — Apt to avoid double counting in the perturbative and 
power correction pieces. Since Apt depends on the coupling this change has an effect 
on the central values of the fit. For example, using the formula of [7, 8, 45] with 
the current data set one obtains at the next-to-leading order (A; = 2) a central value 
of af^{Mz) = 0.124 (the Milan factor [45] is not included) rather than af^ (Mz) = 
0.128. 
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Comparing tabic 3 with the resummation^ results of table 2 we find that the 
truncation leads to an overestimated value of Us- This comparison invalidates the 
next-to-leading order procedure of ref. [7, 8, 45]. As more terms are included, the 
value of as becomes closer to the resummation result, af^ (Mz) = 0.111. This value 
is not reached even for A; = 6: indeed there is a difference between a fixed order 
calculation and a regularized sum, e.g. in the principal value regularization. The 
latter is close, as we saw, to truncation of the series at the minimal term, but then 
the order of truncation is not fixed but rather depends on Q. 

The most drastic change in table 3 is, of course, between the next-to-leading 
order based fit {k = 2) and the /c = 3 fit which includes an estimated (5.21) next-to- 
next-to-leading order contribution from the "leading skeleton", a|is with — 188. 
There is some uncertainty in the estimated coefficient ^2 so long as the identity of 
the "skeleton coupling" is not known. Here we assumed that Cpx is in the "gluon 
bremsstrahlung" scheme, and so we used di = 1 — ^ in eq. (5.21). If we assume 
instead^ the pinch technique scheme, with di — 1, we obtain — 279, which yields 
a best fit for A; = 3 at af^ (Mz) = 0.114 with Apert = 0.66 (cf. A; = 3 in table 3). Note 
that if di (which characterizes the relation between the "skeleton scheme" and MS) 
is not large, it is the "large /5o" term (which is proportional to /9q and independent 
oidi) that dominates in eq. (5.21): t9, = 1.5776 dj + 37.59 + 239.62. 

We stress that the choice we made in (5.20) to expand if^lfi^) in terms of (Q^) 
is arbitrary. We could in principle pick any renormalization scale and scheme. A 
particularly good choice, using still the MS scheme, is to set the scale equal to the 
BLM scale, namely to eliminate the term proportional to Pq from the next-to-leading 
order coefficient 

/xLM = Q'exp(^^ + cij (5.39) 

yielding //g^M — 0.0447(5. We then obtain from (5.20) the following series^ 

= fo ajis (/^blm) + di fo (/^Lm) (5-40) 
, / - 2 . 1 3/i^ + 7rVo^-3/2/o 2I ^3 / 2 N , 

+ I "1 /O - g Po f «MS IA^BLmJ H 

which provides a good approximation to -Rapt — -Rpt|pv already at the leading order, 
as shown in fig. 7. We mention that a particularly low renormalization scale was 
suggested for this observable in [12, 13]. Such a choice can now be justified from 
another view point, noting that the BLM scale approximates well the resummed 
perturbative series. 

^We recall that the coefficients in (5.20) are computed based on the one-loop a^g and so the 
relevant comparison is with the one- loop resummation fit in tabic 2. 

^The dependence of the resummation based fit on the "skeleton scheme" is discussed in sec. 5.6. 

II Since ci can be swallowed into the definition of the scale, it disappears from the BLM series 
(5.40) completely. Note also that in this expansion all the higher order terms linear in /3o vanish. 
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The proximity of the leading order BLM result to the Borel sum in fig. 7 sug- 
gests that performing the fit based on a next-to-leading order partial sum in MS with 
= AiBLM would be much better than with Hr = Q corresponding to A; = 2 in table 3. 
Indeed, performing such a fit (with a power term of the form Xpert/ Q), we find a sig- 
nificant change in the extracted parameters. The central value is af^ (Mz) = 0.116 
(with Apert = 0.55 GeV), which is much closer to the best fit result of our resumma- 
tion, af^ (Mz) = 0.111. It should be noted that the results do not coincide: leading 
order BLM scale-setting is not a substitute to actually performing the resummation 
(see related observations in [25]). 

Note that the success of BLM in MS is not guaranteed a priori. It is the smallness 
of di in the relation (2.7) between the scheme coupling (chosen as MS) and the 
"skeleton coupling" (assumed to be the "gluon bremsstrahlung" coupling) which plays 
a role here. The most natural scheme to apply BLM is the "skeleton scheme" CpT, 
where ci = di = 0. Then, similarly to the Euclidean case (2.15), there is a scheme 
invariant interpretation to the BLM scale 

/^LM = Q'exp(^|j (5.41) 

as the center of ^(e). In the case of the average thrust it is 

//BLM = 0.1028 g, (5.42) 

as can be verified directly in fig. 4. Now the leading term in the BLM expansion 
(5.40) is /o apT(A''BLM) ^^'^ the next-to- leading order correction vanishes. Higher order 
corrections have a simple interpretation [26, 29] in terms of the properties for the 
function T{e). For instance, the next-to- next-to-leading order is related to the width 
of J-'{e) through its second moment f2- 

5.5 Cutoff reguleirization and the infrared finite coupling approach 

As explained in sec. 3.3 the cutoff regularized perturbative sum R^'^ of eq. (3.39) is of 
special interest because it is fully under control in perturbation theory. In particular, 
in this regularization the replacement of the all order coupling Opx by a one-loop or 
two- loop coupling is justified provided /i//A is large enough. In addition, as we saw 
in sec. 4, the infrared cutoff regularization appears naturally in the framework of 
the infrared finite coupling [7, 8] . Let us therefore repeat the analysis of the average 
thrust in terms of this regularization. 

In general, B^y can be obtained from -Rapt by R^^ — Rapt — Ai?. Since in our 
case Re{5i?APT} is negligibly small, we simply have AR ~ Re {Rj^}- The 0{1/Q'^"') 
contribution corresponding to the n-th order term in (5.14) is given by eq. (3.71) 
and eq. (3.98) with (3.85) (or equivalently eq. (3.102)) in the one and two loop cases. 
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respectively. The leading n = ^ term is given in the one-loop case by 

(") 

Rek^a = fel^^e-^Eif^) (5.43) 
with tj = In (yuf/A^) and in the two-loop case by 

^^{^-.i)L, - (t) (4)' ' h-l)} '-^) 

with ij = l/{Podi) where aj is defined in (3.84) and S = 6i = ^ iPi/Po) according to 
eq. (3.83). 

Choosing = 2 GeV, with the same value of a<j as in fig. 5, we obtain the cutoff 
regularized sum (t)p^ = ^R^y which is presented in fig. 9. Comparing fig. 9 with 
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Figure 9: Cutoff regularized perturbative sum with ^u/ = 2 GeV, given a^^ (Mz) = 0.117. 
As in fig. 5, the resummation is calculated with ap^ as a one or two loop coupling in the 
"gluon bremsstrahlung" scheme. In the two-loop case we show also a curve that includes 
the full 0{a1) term, as well as a best fit curve that includes in addition a fitted X/Q 
correction (5.45). 

fig. 5, the first observation is that R^ylQ"^) is significantly lower than -Rapt(Q^) (or 
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-RpT|pv(Q^))) especially** for low Q. This means that a major contribution to the 
resummation (3.2) is from space-like momentum scales below the cutoff fij = 2 GeV, 
i.e. from large distance scales in which the perturbative treatment does not apply. 
This point will be elaborated in the conclusion section. 

As discussed in sec. 3.3, the replacement of the all order coupling by some ap- 
proximation like the one or two loop coupling is best justified in terms of the cutoff 
regularization. It is thus important to verify that i?^v remains stable going from one 
to two loops. Fig. 9 shows that indeed even at low Q the difference between the one 
and two loop R^y{Q'^) is rather small. We conclude that R^y{Q'^) is under pertur- 
bative control for fij = 2 GcV. To improve further the accuracy at which K^^^Q'^) is 
determined one can, in principle, either go to higher orders in the f3 function of dp^ 
or increase /x/. In the latter case, however, one should be careful not to invalidate 
the expansion of T{iJi? /Q"^) in Ai? (or ^(/c^/Q^) in i?") which holds only for small 
enough //f/Q^- The physical value of separating at once short distance physics 
from large distance physics and the perturbative domain from the non-perturbative 
one, should be set satisfying these two constraints. 

Next, we fit the data with 

(0 = ^Kv + 5..o] + ^. (5.45) 

The best fit values of A2GeV for ^ms (Mz) = 0.117 (fig. 9) are listed in table 4 below. 
Note that the fit (5.45) in fig. 9, and the fit (5.35) in fig. 5 are almost identical. This 
is reflected also in the similar values quoted in tables 4 and 1. The reason is that 
the two regularizations, -Rapt (or i?pT|pv) and R^{Q'^) differ, to a very good approx- 
imation, just by the 1/Q power correction of the type we add as a free parameter 
in the fit (see sec. 4). One can thus check that the difference between the fitted 
parameters in the two regularizations is as implied by the calculation. For instance, 
in the one-loop case, the difference is calculable from eq. (5.43), namely 



«7 

eft 


A2 Gov (GeV) 


/point 


one — loop 


1.32 


2.43 


two — loop 


1.14 


3.46 



Table 4: Power term (A2Gcv) fit results with af"" (M^) = 0.117 based on (5.45). 

**Fig. 9 shows that R^{Q'^) is not a monotonous function of Q. This is a unique feature of this 
regularization (with a reasonably high jij) which distinguishes it from both -Rpt|pv(Q^) a.nd typical 
truncated perturbative series, which are monotonic decreasing functions of Q. 
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The sub-leading power term making i?pT|pv and R^y different is related to the 
next non- analytic term in the expansion of J^(e), n = |, which leads to negligibly 
small power corrections. As mentioned above, the latter are calculable using 

eq. (3.71) and (3.102) in the one and two loop cases, respectively. 

The effect of leading and sub-leading power terms in the expansion of ^(e) 
in (5.14) is summarized in table 5 for the lowest experimentally relevant energy 
g = 12 GeV. 



n 


^Re{R--} 


^Re{5i?APT} 




^ (i?- - i?-) 


Cf daft 
2 -"-< 


1 

2 


0.0796 





0.0796 


-0.0401 


0.1197 


1 





0.701 10-3 


0.701 10-=^ 


0.0341 


-0.0334 


3 
2 


-0.131 10-2 





-0.131 10-2 


-0.666 10-2 


0.535 10-2 


2 


-0.822 10"^ 


0.663 10"^ 


-0.821 10-^ 


0.353 10-^ 


-0.117 10-3 



Table 5: Contributions 0{1/Q'^'^^) from increasing order terms in the expansion of J-{e) to 
the difference between different regularizations of the perturbative sum with the one-loop 
running coupling in the "gluon bremsstrahlung" scheme, at Q = 12 GeV. The infrared 
cutoff is set to m = 2 GeV. For comparison, at this energy {t) ~ ^Rapt = 0.1233. 

Table 5 contains contributions of both analytic and non-analytic terms in ^(e). 
It presents the two possible separations of AR discussed in sec. 3.3 - the separation 
into Re{i?fj^} and Re{5i?APT}, corresponding to eq. (3.69), on the left side of the 
table, and the separation into R^^ — R7 and i?^^^, corresponding to eq. (3.67) on 
the right. We find that in the separate pieces on the right the sub-leading terms (e.g. 
n = 1) are relatively important, while not so on the left, where the total AR can be 
approximated by the leading Re{i?"} piece 

Ai?^Re{<J. (5.47) 

Next, consider the dependence of the cutoff regulated perturbative sum Rly{Q'^) 
on the cutoff scale Any two cutoff regularizations differ, in general, by renormalon 
related infrared power corrections. In our case, they differ to a very good approxi- 
mation by a 1/(5 term which can be calculated using eq. (5.43) and (5.44) in the one 
and two loop cases, respectively. As a result, one would practically ol:)tain the same 
best fit in (5.45) for any arbitrary A < ni <^ Q. In the two-loop case, R^^'^{Q'^) is well 
defined even for the extreme choice fij = A. The resulting regularization is shown 
in fig. 8 together with the standard ni = 2 GeV choice. The /x/ = A regularization 
turns out to be close to the best fit in this case. This regularization is, of course, 
unacceptable in the infrared finite coupling approach [7, 8], where one requires that 
the non-perturbative coupling a is well approximated by above iij. As explained 
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in sec. 4, within the infrared finite couphng approach, the non-perturbative param- 
eter A^^ acquires a physical interpretation as the small gluon virtuality moment of 
the "skeleton coupling". Using (4.14) with (5.14), we have 



We stress that the identification of the fit parameter A^^^ with a small virtuality mo- 
ment of the coupling is meaningful only under the strong assumption of universality: 
it is the same "skeleton coupling" for any observable. Eq. (5.49) as written implies 
that the coefficient in front of the "skeleton coupling" average is simply related to the 
small /x^ behavior of the characteristic function, in spite of the non-inclusive nature 
of the thrust. In ref. [45] it is argued that this coefficient is in fact modified by the so 
called "Milan factor" arising at the next-to-leading order level, taking into account 
the emission of two gluons. We ignore these corrections here (see sec. 5.6). 

Let us now extract A2GeV based on the best fit of (5.45). As we know al- 
ready from sec. 5.3 (see table 2), the best fit is obtained in the one-loop case at 
af^ (Mz) = 0.111. The corresponding power correction coefficient in the cutoff regu- 
larization is A2GeV — 1-57 ± 0.12 GeV, where the error corresponds to the combined 
systematic and statistical experimental uncertainties. In the two-loop case, the best 
fit is at af^ (Mz) = 0.110 with AsceV = 1-49 ± 0.12 GeV. In the latter case the cutoff 
regularized sum and the best fit are shown in fig. 10. 

According to eq. (5.49) we have for /jLj — 2 GeV 



which implies 




(5.49) 




one — loop 
two — loop 



(5.50) 



where a is the "gluon bremsstrahlung" coupling. 
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Figure 10: Cutoff regularized perturbative sum with ^/ = 2 GeV, for the best fit value 
a^^ (Mz) = 0.110. As in fig. 6, the resummation is calculated with as a one or two 
loop coupling in the "gluon bremsstrahlung" scheme. In the two-loop case we show also a 
curve that includes the full 0{a'^) term, as well as a best fit curve that includes in addition 
a fitted X/Q correction (5.45). 

5.6 Sources of theoretical uncertainty 

The experimental error on the value of and A was discussed in sec. 5.3. The 
theoretical error, which is harder to determine, will be briefly discussed here. 

A primary source of uncertainty is the unknown identity of the running coupling 
used in the calculation of -Rapt- In sec. 5.3 we approximated this coupling by the 
one or two loop coupling in the "gluon bremsstrahlung" scheme with j3o and (3i 
calculated for Nf = 5. In principle, this coupling should be an "all-order" coupling 
in the "skeleton scheme" (yet unknown), where the number of active flavors depends 
on the scale. Let us try to quantify the errors. 

Five active flavors: Considering the center of mass energy of most relevant exper- 
iments, using* the running coupling with Nj = 5 is known to be a reasonable 

*Another source of error is the assumption, taken in the calculation of J-' and the perturbative 
coefficients (5.3), that the primary quark and anti-quark are massless. This problem is common 
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approximation in the standard pcrturbative approach. However, in the re- 
summation performed here the couphng runs over small gluon virtualities as 
well, and then there are only four or three active flavors. The significance 
of the Nf = 5 approximation in our approach can be checked by calculating 
the difference between the contribution to perturbative sum it^APT from small 
space-like momentum scales with A^^ = 3 or 4 and the corresponding contri- 
bution with Nf — 5. We use for this calculation'^ the formula for R^^J at the 

J IR,2 

one-loop order, which is obtained in the Appendix, eq. (A.4) and (A. 5). The 
relative difference between R^^l_{Nf = 3) and R^^1{Nf = 5) for the lowest 
experimental data point Q = 12 GcV, where the effect is the largest, is just 1.9 
percent. This is negligible compared to other sources of error. 

Two-loop coupling: In the phenomenological analysis we used the one or two loop 
couphng to evaluate the perturbative sum, instead of the "all order skeleton 
coupling". As explained in sec. 3.3 this replacement is justified in the cutoff 
regularization provided /i//A is large enough. Indeed wc saw in sec. 5.5 (see 
figs. 9 and 10) that the cutoff regularized sum with /i/ = 2 GcV does 
not change much going from one to two loops. This stability suggests that 
the replacement of the all order coupling by some low order coupling is a 
reasonable approximation for this value of the cutoff. To go further one would 
like to examine the effect of a three-loop correction P2 in the P function of the 
"skeleton couphng" on R^'^. However, technically, this is a complicated task. 
Below, we shall examine instead the effect on i?APT- At the one and two loop 
level i?APT is rather stable, as shown in figs. 5 and 6. A priori, the stability 
of -Rapt or the principal value Borel sum may be less intuitive because these 
regularizations do involve the coupling in the infrared. The stability of -Rapt 
may be explained in terms of the infrared stability of the APT coupling itself, 
which is discussed in [32, 40, 38], and goes beyond the two-loop level. Returning 
to our case, the effect of a three- loop correction in (3 (a) on i?APT can be examined 
explicitly. Since we do not know the identity of the "skeleton coupling" , we just 
try various choices of P2 where for simplicity we use the Pade improved form of 
the three-loop (3 function which can be solved analytically using the Lambert W 
function [38]. As an example consider the case where the f3 function is modified 
to include a large three-loop term, f32/(3o = —45 (both positive and negative 
values were studied). Taking the world average value of {af^ (M^-) = 0.117) 
we obtain at the lowest experimental data point Q = 12 GeV, where the effect 
is the largest, the following results for {t) ~ ^Rapt- at one-loop {t) ~ 0.123, at 
two-loop (t) ~ 0.131, and at three-loop (t) ~ 0.116. These differences are small 

to the standard perturbative approach and to ours and it certainly deserves attention. We do not 
attempt to estimate this error here. 

'^One could evaluate also Re < R^Ji \ or ARi which are all roughly the same in this case. 



55 



compared to the experimental error bars. The conclusion remains the same: 
the perturbative sum is quite stable with respect to including higher orders in 
the /5 function. Whichever regularization is chosen, what counts at the end is to 
what extent the extracted parameters depend on the approximation used for the 
perturbative coupling. It is clear from table 2 that the resummation procedure 
is quite stable going from one to two loops concerning the determination of ag- 
In particular, the resulting uncertainty in the extracted value of ag is smaller 
than the experimental error (5.36). Moreover, as we saw, one can use either 
-Rapt oi' -^uv obtaining the same value of Ug- Thus, based on the stability 

of -^APT 

at three-loop we conclude that the extracted value of will remain 
stable, despite our lack of knowledge concerning the stability of at this 
level. On the other hand, for the determination of the power correction or the 
low momentum average of the coupling, the relevant theoretical uncertainty 
is comparable to the experimental error. This is reflected for example in the 
difference between the one and two loop results in tables 1 and 4. Similarly we 
have for* ag — 0.110, 





APT / PV 
A (GeV) 


cutoff 
A2 GeV (GeV) 


one — loop 


0.79 


1.62 


two — loop 


0.62 


1.49 



Table 6: Power term fit results with af^ (M^) = 0.110. 

Here it is important, in principle, which regularization is used. In our case, 
as shown in table 6, the cutoff regularization with /ii — 2 GeV yields slightly 
more stable values of the power correction coefficient than the principal value 
Borel sum (or APT) regularization. Nevertheless, this uncertainty in A^^ is 
significant and should be taken into account when discussing the universality 
of the infrared finite coupling. 

"gluon bremsstrahlung" scheme: In the phenomenological analysis above we 
used the "gluon bremsstrahlung" renormalization scheme as the scheme of the 
"skeleton coupling" apT. As explained in sec. 2, the "skeleton coupling" should 
actually be defined diagrammatically in a unique way and the only thing we 
know about it, is that it coincides with the V-scheme in the Abelian limit. 
Alternative schemes were discussed in sec. 2 and 5.2. Of particular interest is 
the pinch technique coupling which has some systematic justification [22] and 

■l-Note that this value of as is close but not equal to the best fit value in the one-loop case, 
table 2. 
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which differs from the "gluon bremsstrahlung" scheme by having a relatively 
large next-to-leading order correction (5.28) on top of the terms included in the 
resummation. The leading 0{a^) "scheme dependence" effects corresponding 
to taking Opx in different schemes cancel between -Rapt and remaining next-to- 
leading order correction in eq. (5.34). The residual dependence is thus formally 
of order a^. In order to check the sensitivity of our procedure to the choice of 
OpT, we repeated the whole analysis with cip-p as G Starting as in sec. 5.3 
with the world average value of as, af^ (Mz) = 0.117 which corresponds to 
Ai-|^°°p = 0.274 GeV and A^:|X = 0.686 GeV, we obtain the results shown in 
fig. 11. Here the resummation result, namely -Rapt or principal value Borel sum 




Figure 11: Principal value Borel sum (or -Rapt) regularization of the perturbative sum 
(3.2) with the "skeleton coupling" Opx in the pinch technique scheme [22], given a^^ (Mz) = 
0.117. For more details, see the caption in fig. 5. 

with ttpT at one-loop or two-loops, overshoots most data points and especially 
the large Q ones. However, after including the next-to-leading order correction 

(5.34) with the appropriate negative coefficient (5.28), the perturbative result 
becomes again quite close to the large Q data points. Adding a power correction 

(5.35) the best fit (with Op^ at two-loop) is obtained for A = 0.44 GeV, with 
X^/point = 4.35. In the infrared cutoff regularization with fij = 2 GeV the same 
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fit is obtained with A2GeV = 1.53 GeV. Finally, we fit also using Opx in the 
pinch scheme. The results of the best fit are shown in fig. 12. The best fit with 




Figure 12: Principal value Borel sum (or -Rapt) regularization of the perturbative sum 
(3.2) with the "skeleton coupling" in the pinch technique scheme [22], for the best fit 
value of (Mz) = 0.109. For more details, see the caption in fig. 5. 



the two-loop running coupling in this scheme is obtained for af^ (Mz) = 0.109, 
with a power correction A = 0.83 GeV in the principal value regularization and 
A2GcV = 1-85 GeV in the cutoff regularization. In this fit x^/point = 1.37. The 
comparison of these results with eqs. (5.36) and (5.37) suggests that precise 
identification of the "skeleton coupling" Opx is not crucial in this case for the 
determination of a^. 

Another possible source of uncertainty are 0{a^) terms which are not related to 
the "leading skeleton" . At the single dressed gluon approximation level, such effects 
can be taken into account by calculating 

{t),AQ') ^ ^ [sriQ') + iia' + ka'] (5.51) 

instead of (5.17) where only the next-to-leading term was included. The coefficient t2 
can be determined from the full next-to-next-to-leading order perturbative coefficient 
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t2, t2 = t2 — once the time-like "skeleton coupling" af^ in (3.2) and the scheme 
coupling a are specified. Unfortunately, t2 is not known yet. Since the running 
coupling effects we resum here are understood to be the major source of growth of 
the standard coefficients tj, the general expectation is that the remaining coefficients 
ii are small. Nevertheless, let us see what are the consequences if this assumption 
does not hold. We proceed as before using the coupling a — -Rapt(Q^)//o 



I + — a + — a 
Jo Jo 



(5.52) 



For the perturbative expansion to make sense one should require that the next-to- 
next-to-leading order term would be smaller than the leading term, i.e. t2//o < V^^ 
or i2//o < fo/R-APT- For the lowest data points Q — 12 GeV the r.h.s. equals roughly 
64 so we shall consider possible next-to-next-to- leading corrections with t2//o ^ 64. 
Once t2 is fixed we can fit the thrust data with (5.52) plus a X/Q power term. 
The results of such a fit turn out to depend significantly on the value of ^2- For 
instance, with t2//o = —63, the best fit is obtained for af^ (Mz) = 0.117 with A = 
2.4 GeV and xVpoint = 1.478 while for t2//o = 56 the best fit is af^(Mz) = 0.104 
with A = 0.96 GeV and xVpoint = 1.313. Note that one could gain nothing by 
performing a three parameter fit, varying t2 in addition to ag and A, since the 
curve becomes almost flat. The conclusion is that existence of a large next-to-next- 
to-leading coefficient on top of the resummation performed would indeed modify the 
fit results significantly. 

Finally, we relied heavily on the assumption that the inclusive resummation pro- 
gram based on the "skeleton expansion" applies to the case of the average thrust. The 
physical intuition behind this assumption is that most gluon splittings are roughly 
collinear, and thus one can ignore the cases where the split is to opposite hemi- 
spheres. This was further justified in sec. 5.2 by showing that the magnitude of the 
non-inclusive correction at the next-to-leading order is small, at least as far as the 
"Abehan" (/?o dependent) part is concerned. On these grounds we assumed that the 
non-inclusive effect is small also at higher orders. In the large Nf calculation [5] it 
was found that non-inclusive corrections do affect the asymptotic factorial growth of 
the perturbative coefficients, and therefore also the 1/Q power correction. However, 
if such corrections are small before the asymptotic regime it reached (like they are 
at the next-to-leading order) they would not change the extracted value of a^. Fur- 
ther, concentrating on the single gluon emission approximation, we did not take into 
account perturbative and non-perturbative effects which are related to emission of 
more than one gluon. Such effects were analyzed in [45] at the next-to-leading order 
level, and were argued to give a significant contribution to the 1/Q power term. To 
understand how to incorporate this type of corrections into the present approach 
requires further work. 
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6. Conclusions 



In this paper we suggest a calculation method for inclusive enough observables in 
QCD, which is based on the conjecture that the Abelian skeleton expansion can be 
extended to the non- Abelian case. We did not make here an attempt to establish 
the existence of such an expansion in QCD, but rather took a practical approach to 
use it at the single dressed gluon level based on the analogy with the Abelian case. 
Clearly, to go beyond this approximation, a more rigorous treatment is required, 
which would hopefully establish a diagrammatic interpretation to the sub-leading 
terms in the "skeleton expansion" and an all order identification of the "skeleton 
coupling". Such a program was initiated by Watson in [22] based on the pinch 
technique. It certainly deserves more attention. 

The assumption of a "skeleton expansion" leads to the following crucial conse- 
quences already at the single dressed gluon approximation: 

a) resummation of all order perturbative corrections that are related to the run- 
ning coupling (renormalons) and parametrization of power corrections must be 
performed together. 

b) the resummation yields a renormalization scheme invariant result. 

The idea that a certain type of infrared power corrections can be related to a 
universal infrared finite coupling [7, 8] fits very well into the "skeleton expansion" 
approach. If the "skeleton coupling" is well defined on the non-perturbative level 
down to the infrared, there is no renormalon ambiguity and the perturbative sum plus 
the power corrections are obtained at once by performing the renormalon integral. 

In this work we considered the specific example of the average thrust in e'^e" 
annihilation, which is just appropriate to demonstrate the advantages of the sug- 
gested procedure over the standard perturbative treatment. For this observable the 
available next-to-leading order perturbative calculation with a standard choice of 
renormalization scheme and scale turns out to be far from the experimental data 
(fig. 1). The perturbative series suffers from large renormalization scale dependence. 
We assume that the most important higher order corrections are related to running 
coupling effects which have a precise meaning and can be resummed in the "dressed 
skeleton expansion" framework. Assuming a "skeleton expansion" with a uniqely 
identified "skeleton coupling" , once the resummation is performed there is no more 
scale or scheme dependence. The resummation itself, however, is ambiguous due to 
infrared renormalons. The ambiguity is resolved only when the regularized perturba- 
tive sum is combined with the appropriate power correction. Since for this observable 
the leading renormalon is so pronounced (n = |) one cannot ignore the resummation 
ambiguity and the related 1/Q power corrections. 

At the first stage we perform the "leading skeleton" resummation in a standard 
regularization such as the principal value Borel sum (sec. 5.3) or truncation of the 
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series at the minimal term (sec. 5.4). We find (figs. 5 and 6) that a major part of the 
discrepancy between the next-to-leading order result and experimental data can be 
explained in terms of the resummation. The most important consequence is that the 
extracted value of is different than in a standard next-to-leading order based fit. 
We obtain af^ (Mz) = O.llOi 0.002, where the error reflects the combined statistical 
and systematic experimental uncertainties. The discrepancy between this result and 
the world average value of ag is interesting and calls for further study. It is necessary 
to apply a similar resummation program to other QCD observables that are prone 
to having large higher order corrections of the same type, e.g. other average event 
shapes^, in order to extract a reliable value for ag from experimental data. 

At the next stage (sec. 5.5) we specify to the infrared finite coupling approach 
and reinterpret the resummation plus power corrections as originating together in 
the renormalon integral performed down to zero momentum. The regularization we 
use in this case is an infrared cutoff Hi — 2 GeV on the space-hke momentum. We 
assume that the full "non-perturbative" coupling coincides with the perturbative one 
above this scale. Thus the cutoff separates short and long distance physics as well 
as perturbative and "non-perturbative" physics. Using the cutoff regularization, the 
perturbative sum (i?uv) includes only short distance contributions related to scales 
where the "skeleton coupling" is controlled by the leading terms in the perturbative 
j3 function. The "non-perturbative" long distance contribution is expressed in terms 
of moments of the assumed universal "skeleton coupling" a{k'^) in the infrared re- 
gion: < A;^ < /Xj. Comparing fig. 9 and 10 with fig. 5 and 6, we flnd that the gap 
between the cutoff regularized sum and the standard regularizations discussed 
previously (e.g. principal value Borel sum) is quite large, especially for low Q values 
(see also flg. 8). This means that an important part of the contribution to what is 
called perturbative sum^ in the standard regularization actually originates in small 
momentum scales where the coupling is non-perturbative. On the other hand, this 
does not imply that the resummation is entirely related to long distance contribu- 
tion: at large Q, becomes quite close to the standard principal value Borel sum 
regularization. 

It is useful to examine the physical scales which contribute to the perturbative 
sum by applying the ELM scale-setting method to the "leading skeleton" term. For 
Euclidean quantities, the momentum scales from which the major contribution to 
the perturbative sum (say, in the principal value regularization) originates, is the 
BLM scale in the "skeleton scheme" (2.15). This characteristic scale corresponds 
to the average virtuality of the exchanged gluon. It is natural to compare it with 

^The analysis of event shape distributions is more involved since there one has to perform 
multiple soft gluon resummation in addition to the renormalon resummation performed here. 

^Note that the prescription dependence of the separate pieces, namely the regularized pertur- 
bative sum and the power correction, does not appear in the standard treatment when a power 
correction is added to the truncated perturbative sum, e.g. at next-to-leading order. 
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a typical hadronic scale, e.g. in our context with the cutoff /x/. If /Ublm ^ fJ'i the 
resummation program is purely perturbative in the sense that it does not depend on 
the coupling at infrared scales. If /Xblm ^ l^i, the most important contribution to the 
"perturbative sum" comes from long distances and then the role of power corrections 
will be crucial. Applying the same considerations to the Minkowskian representation 
(3.2) we find that it is the center of ^(//^/Q^) which determines the BLM scale in 
the "skeleton scheme" (5.41), the scale where the most significant contribution to 
the perturbative sum comes from". Comparing /Xblm of (5.41) with the (space-like) 
cutoff Hi we can learn to what extent the corrections are controlled by perturbation 
theory. In the case of thrust, ^Iblm in the "skeleton scheme" (5.42) is 0.1028(5, which 
is above the fij = 2 GeV cutoff for most of the relevant experimental data. For the 
low Q data points, the proximity of the two scales implies that an important part 
of the contribution comes from the non-perturbative large distance regime, as we 
learned already based on fig. 9 and 10. 

A crucial assumption we made in this work is that the inclusive "skeleton expan- 
sion" approach is useful also for not-completely-inclusive quantities like the average 
thrust. As explained in sec. 5.1, the non-inclusive nature of the thrust appears first 
in perturbation theory due to a split of a gluon into the two opposite hemispheres 
defined by the thrust axis. This possibility is ignored in the resummation we per- 
form: non-inclusive corrections do not fit into a "skeleton expansion". Physically, 
the assumption that the non-inclusive corrections are small seems reasonable since in 
QCD one expects most parton splittings to be roughly collinear. We further justified 
it in sec. 5.2 by showing that the "Abehan" (/5o dependent) part of the next-to- 
leading order term that emerges from the inclusive gluon dressing almost coincides 
with the corresponding term in the full next-to-leading QCD calculation - the two 
would coincide exactly for an inclusive quantity. It would be interesting to examine 
how well this approximation works at higher order, i.e. what is the error due to ne- 
glecting the "Non — Skeleton" terms in eq. (5.16) beyond the next-to-leading order. 
The Abelian case can be studied to all orders in perturbation theory comparing the 
present inclusive "massive gluon" calculation with that of the non-inclusive large Nf 
renormalon calculation [5] , similarly to the comparison made in [48] for the longitu- 
dinal cross section. Further contributions relevant to the 1/Q power term which are 
specific to the non- Abelian theory were analyzed in ref. [45] at the next-to-leading 
order level. The connection of this type of corrections with the "skeleton expansion" 
and renormalons should be clarified. 

It is tempting to speculate about the perturbative nature of the "non pertur- 
bative" infrared "skeleton coupling" , namely that the infrared part of the "skeleton 
coupling" could itself be determined from information contained in perturbation the- 

I'We note, however, that the integral (3.2) runs over time-like momentum, and so the BLM scale 
in (5.41) cannot be identified immediately as related to some characteristic distance scale. 
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ory. It could be, for instance, that the all order Borel resummed "skeleton coupling" 
is infrared regular, i.e. free of Landau singularities**. This last possibility seems 
natural as a continuation of the perturbative fixed point scenario which presumably 
holds within the conformal window [49]: for 10 < Nf < 16 there is a fixed point in 
the perturbative /3 function and the perturbative coupling appears to be causal. One 
might then speculate that analytic continuation in Nf below the conformal window^^ 
might yield the correct infrared finite "skeleton couphng" . In this sense the physics 
of the universal infrared finite "skeleton coupling" , and in particular that of the gluon 
condensate, would be entirely perturbative. 

Finally we recall that the BLM scale-setting procedure is aimed to approximate 
the "leading skeleton" integral. In the case of the average thrust we saw that choosing 
the BLM scale leads to a significant improvement of the next-to-leading order partial 
sum compared to the naive choice of scale /ir — Q- However, a fit based on the 
next-to-leading order series with hr — //blm, gives af^{Mz) — 0.116, which is still 
appreciably different from the resummation based fit. 
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**For a realistic number of flavors, the first two universal coefiicients of the /3 function are negative. 
It is therefore quite clear that without resummation a perturbative fixed point cannot appear in 
any coupling, including the "skeleton coupling". Sec the related discussion in [37, 38, 49]. 

^^This does not imply that more genuine non-perturbative phenomena like spontaneous chiral 
symmetry breaking could be studied by analytic continuation in Nf from the conformal window, 
since there these phenomena are probably absent. 
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A. Cutoff regularization in terms of APT coupling 

One can give a useful alternative expression for AR, based on the general formula 
eq. (4.10) in sec. 4 specialized to the APT couphng 



AR = + SRf^-" 



where each of the terms is separately unambiguous*-!- . 
The first term 

can be evaluated exphcitly in the one- loop case (3.14) 



(A.l) 



(A.2) 
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Consider a generic term in the small e expansion of J-{e), as in (3.27). Similarly to 
eq. (3.48) we have 
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In the one-loop case these integrals give 

^nil-loop = - e""*^ Eil {-nti) - Lerch$ (e*^ 1, n) 







and 



(A.4) 

(A.5) 
(A.6) 

(A.7) 



ntie-''*' Eil {-nti) + 1 - nLerch$(e*^ 2, n)! (A.8) 



-loop 

where as before tj = log(//f/A^), Eij^ is defined in (3.55) and Lerch$ is defined by 
analytic continuation of 



Lcrch$(z, fc, n) = 



(A.9) 



ttEq. (A.l) actually defines SR^^'^. 
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from |2;| < 1 to the whole complex z plane. We have used the property that for a 
positive n, 



z 



n 



dz 



1 — Z Z 



= ZjLerch^{zi, l,ri) 



(A.IO) 



(here zj — A*f/A^) and taking the derivative with respect to n, 

f^i , , z'^ dz 
/ ^og{z)- 

^0 



z z 



log(2;7) Zj Lerch$(2;/, 1, n) - z^ Lerch$(2;/, 2, n). (A.ll) 



In (A. 7) and (A. 8), as opposed to /„ and in (3.53) and (3.54), there is no cut for 
Hi > A since the imaginary part from the Ei function is cancelled by the one from 
the Lerch$ function. This is easy to understand, since the Lerch$ part originates 
in the integration over ^Oapt of eq. (3.14). The simple pole in ^Oapt cancels by 
construction the simple Landau pole in the perturbative one-loop coupling, and thus 
also the related ambiguity in /„ and should cancel. 

One can write explicit formulae for the integral of eq. (A.IO) in terms of elemen- 
tary functions for simple cases such as integer n values 



Lerch$(2;/, 1, n) 



zr'^ ln(l - zi) 



or n 



Lerch$(^;7, 1, ^) = ^ 

2 Jzi 



In 



1 - 



(A.12) 



(A.13) 



1 + 

An interesting empirical observation is that J„ ~ Re{/„} and J\ ~ Re{/^}. The 
ratios Jre/Re{/„} and J^/Re{/^} approach 1 os jii/h. increases, due to the smallness 
of the Lerch$ piece in J„ and J^. Already* at iij/A^ 10 these ratios are not too far 
from 1, as shown in table 7. 



n 


lii/A 


J„/Re{4} 




1 

2 


5 


1.11 


6.33 


1 
2 


8 


1.05 


2.60 


1 


5 


0.72 


0.84 


1 


8 


0.81 


0.80 


3 
2 


5 


0.73 


0.64 


3 
2 


8 


0.85 


0.75 


2 


5 


0.77 


0.66 


2 


8 


0.88 


0.80 



Table 7: Comparison botwven ./„ and Re{/„} and between and Re{/^} for /x//A < 10. 



*An exception is the log term J^/Re {In} for n = ^. 
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This implies that even for not too large /U//A the contribution of a generic non- 
analytic term in the small e expansion of J-'r admits 

Re{<4^i?f-. (A.14) 

If we assume in addition that (3.75) holds, i.e. that Re {5/2^^"^} is negligible compared 
to Re {-Rfij^n} then in eq. (3.69) Ai?„ ~ Re {R^^^J^, which imphes 

^ it:f,^^. (A.15) 

By comparison with eq. (A.l) it follows that already for not too large /i//A, the 
corresponding ultraviolet contribution in 6R^y\ is small. For analytic terms, the 
infrared contributions Re |-Rfj^„| and Rf^^ vanishes identically. 

Note that eq. (A.15) is equivalent to the statement that the /ii cutoffs in the 
Euclidean and in the Minkowskian representations are practically equivalent for the 
APT couphng, in the sense that 

-^nCn ~ ^iR,n — ^<',n ~ ^^<,n (A. 16) 

One can also compute the Minkowskian "ultraviolet correction" 5R^^^ in terms 
of the associated Euclidean one 5D^^ using the Minkowskian-Euclidean connection 
pointed out in [16]. We have 

with 

/•OC flip' 

5D^l^{Q^) = j^^ — 6a^Mk') MkVQ') (A. 18) 

This observation allows to express SR^^'^ in terms of integrals over ^ciapt- Indeed 
one can derive [16] from eq. (A. 18) the large behavior of 6D^^'^{Q'^). Assuming 
e.g. the leading term in the small gluon mass expansion of JF^ is a n = | pure power 
term (this example is relevant to the thrust case), one finds 

SDZ-{Q') ^ -lci^^K{^^J)^ + ... (A.19) 

where 

KM ^ ^ S-a^Mk') (^^j (A.20) 

is ultraviolet convergent, since (JoaptI^^) is 0{l/k'^) at large k^, and may eventually 
be computed analytically for the one and two loop couplings. Eq. (A. 17) then gives 

SRtViQ') ^ -Ici^^Ki^ij) + (A.21) 
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